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ABSTRACT 

A  model  trlton  has  been  studied.   The  model  consists 
of  three  bodies  interacting  via  a  two-body  central  gausslan 
potential.   A  type  of  polarized  gaussian  wave  function  was 
used  in  a  minimal  principle  for  the  energy.   This  polarized 
function  gives  significantly  lower  energy  than  the  unpolar- 
ized  wavefunction.   It  also  has  the  advantages  of  simplicity 
and  generality.   The  binding  energy  compares  very  well  with 
an  exact  numerical  solution  of  the  problem,  as  well  as  the 
energy  obtained  by  other  means.   Similar  results  were 
obtained  for  the  alpha  particle.   A  three  body  potential 
can  be  written  down  in  various  ways  so  that  the  wave 
function  is  an  exact  solution  of  the  Hamiltonian  formed 
by  adding  this  potential  to  the  original  Hamiltonian. 
One  way  of  considering  this  function  is  in  coordinate  space. 
A  simple  function  has  been  found  to  approximate  this  in 
regions  where  the  wavefunction  is  not  small.   It  has  the 
form  of  a  gaussian  times  a  quadratic  expression.   Another 
way  is  by  using  projection  operators.   This  has  been  used 
for  the  following  calculations.   An  upper  bound  for  the 
mean  square  radius  was  calculated.   This  involved  an 
extension  of  the  ideas  of  a  previous  calculation  of  an 
upper  bound  for  the  kinetic  energy  of  two  electrons.   A 
variational  principle  for  quantities  other  than  the  energy 
was  set  up.  This  theory  was  applied  to  the  mean  square  radius 
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VARIATIONAL  STUDIES  OP  THE  FEW  BODY  NUCLEUS 
S.  Aranoff  and  J.  Percus 

1 .   Introduction. 

A  system  of  two  particles  Interacting  via  a  central 
potential   can  easily  be  solved  for  most  potentials,  as 
the  problem  can  be  reduced  to  an  equivalent  one-dimensional 
problem.   The  only  non-trivial  coordinate  is  the  scalar 
distance  between  the  particles.   The  three  body  problem, 
i.e.,  a  system  of  three  particles  interacting  via  a  central, 
two-body  potential,  usually  cannot  be  solved  classically, 
let  alone  quantum  mechanically.   This  is  because  there  are 
three  non-trivial  coordinates  for  the  problem  --  the  scalar 
distance  between  two  particles,  the  scalar  distance  between 
the  third  particle  and  the  center  of  mass  of  the  first  two, 
and  the  angle  positioning  the  third  particle  relative  to 
the  first  two.   The  procedure  followed  in  the  past  for  the 
triton  has  been  to  approximate  the  spatial  part  of  the 
wave  function  in  ways  which  essentially  reduce  the  problem 
involving  three  non-trivial  coordinates  to  a  problem 
involving  only  one  non-trivial  coordinate.   One  way  is  to 
approximate  the  correct  wave  function  by   functionswhich 
area  product  of  two-body  functions,  i.e.,  functions  which 
depend  upon  the  distance  between  two  particles  only. 
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Peshback  and  Rublnow  and  MacDonald  have  studied  the  best 

types  of  wave  functions  which  are  products  of  two  body  terms. 
These  functions,  however,  are  not  forms  which  are  easy  to 
handle  analytically.   Another  way  Is  to  use  as  approxima- 
tions functions  which  depend  upon  quantities  invariant 
under  particle  exchange,  such  as  the  sum  of  the  squares 
of  the  particle  distances.   A  function  which  is  usually 
chosen  is  a  gausslan,  for  then  most  integrals  can  be  done 
exactly.   A  simple  gausslan,  however,  gives  a  very  poor 
result  for  the  binding  energy. 

In  this  paper,  a  procedure  is  set  up  for  approximat- 
ing the  three  body  problem  in  a  fashion  permitting  some  of 
the  distinctive  features  of  the  three-body  problem  to  be 
seen.   The  approximate  wave  function  is  written  as  a  sum 
of  wave  functions  each  of  which  Involves  the  position  of 
the  1-th  particle  in  an  asymmetrical  fashion.   The  procedure 
is  applied  to  the  gausslan,  to  get  an  improved  trial  wave 
function  for  a  calculation  of  the  energies  of  the  trlton 
and  alpha  particle.   The  concepts  of  polarization  and 
resonance  are  utilized  in  this  wave  function, 
permitting  an  easy  generalization  to  other  problems. 
Various  techniques  have  been  set  up  and  applied  to  the 
three  body  nucleus  using  this  polarized  gausslan  wave 
function.   The  mean  square  radius  was  studied  along  with 
the  energy.   The  results  are  very  good,  suggesting  that 
this  polarized  wave  function  is  an  excellent  approximation 
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to  the  correct  wave  function,  and  so  could  well  be  used 
in  studies  of  the  few  body  nucleus.   For  phase  shift 
calculations.  In  particular,  not  only  Is  such  a  function 
simple  to  handle,  but  has  the  additional  advantage  that 
the  distortion  of  the  scatterer  during  the  scattering 
can  easily  be  taken  into  account. 

In  Section  2,  the  model  of  the  nucleus  is  set  up. 
In  the  model  trlton  or  helium  nucleus  studied  below, 
three  or  four  nucleons  Interact  via  a  central,  two-body 
gausslan  potential,  with  all  other  factors  such  as  hard 
core,  tensor  and  spin  forces,  and  relativlstlc  effects 
neglected.   These  approximations  allow  us  to  study  the 
essential  ideas  without  undue  complications.   In  Section  3^ 
the  polarized  wave  functions  which  minimize  the  energies 
of  the  trlton  and  alpha  particles  are  exhibited  and  discussed 
In  Section  4  various  forms  of  a  three-body  correction 
potential  are  discussed.   This  potential  is  useful  for 
studies  of  the  difference  between  the  approximate  wave 
function,  i.e.,  the  polarized  wave  function,  and  the  true 
solution  of  the  problem.   This  correction  potential  is 
defined  by  the  property  that  the  polarized  wave  function 
is  an  exact  eigenfunctlon  of  the  hamiltonlan  formed  by 
adding  this  correction  potential  to  the  original  hamiltonlan. 
In  Section  5,  a  method  for  finding  an  upper  bound  on  the 
error  in  the  mean  square  radius  is  discussed.   In  Section  6 
a  variational  principle  for  quantities  other  than  the  energy 


Is  set  up .   The  theory  is  then  applied  to  the  mean  square 
radius,  with  excellent  results. 

2.   The  Model  Studied. 

We  consider  a  nuclear  model  consisting  of  three  or 
four  nucleons  with  a  two-body  gausslan  potential,  with 
no  tensor  or  spin  forces  and  no  hard  core.   Coulomb 
forces  are  considered  for  the  helium  nuclei.   Since  no 
spin  forces  are  considered,  all  spins  occur  with  equal 
probability,  the  experimental  fact  that  the  deuteron 
has  spin  1  notwithstanding.   An  additional  consequence 
of  ignoring  spin  forces  is  that  the  total  spin  s  and 
the  fotal  angular  momentum  i  are  good  quantum  numbers. 
We  therefore  have  the  lowest  state  with  jl   =  0 .      This  is 
possible  as  the  total  wave  function  must  be  totally 
antisymmetric,  and  the  product  of  spin  and  isotopic  spin 
for  four  or  less  nucleons  can  be  antlsymmetr-ized  with 
the  space  part  symmetric  with  respect  to  permutations 
of  the  nucleons .   It  is  shown  by  the  following  argument 
that  the  problem  reduces  to  a  hamiltonian  operating  upon 
a  space  symmetric  wavefunction.   The  Interaction  between 
two  nucleons  has  been  written  as 

(/(r.  .)  =  [wW.  .+mM.  .+bB.  .H-hH.  .]V(r.  .) 

3  4 
by  Laskar^'  ,  where  W.  .  is  the  identity  operator,  M.  ., 

B.  .,  and  H.  .  are  the  Majorana,  Bartlett,  and  Heisenberg 

operators  exchanging  space,  spin,  and  spin-and-gpace 
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coordinates.   In  taking  expectation  values  of  U{v.  .) ,    the 
terms  b  and  h  disappear  as  the  wave  function  contains  both 

the  symmetric  and  antisymmetric  spin  functions  equally,  as 

2 

.—^    ~ur  •  • 

there  are  no  spin  forces.   We  set  V  =  Vq  >   e  '^^ ,    and 
normalize  w+m+b+h  =1,  so  that  V^  Is  the  depth  of  the 
well  of  the  deuteron  In  the  -^S  state,   x  =  w+m-b-h  Is 
then  the  ratio  of  the  S  to  the  ^S  Interaction.   Laskar 
has  taken  x  =  0.6.   We  get  w+m  =  -^(1+x)  =0.8.   For  a  space 
symmetric  function,  M.  .  Is  the  Identity  operator,  and  so  we 
get'Mr.j)  =  .8v(r.j). 

Laskar  has  chosen  the  potential  to  be  what  he  felt 
to  be  the  best  effective  central  potential  for  the  light 
nuclei  (two,  three,  and  four  body  nucleus),  which  replaces 
the  actual  complicated  forces.   He  has  arrived  at  this  by 
following  Bransden-^   to  get  the  range  \i ,    and  by  Interpolat- 
ing the  work  of  Burke  and  Robertson   to  get  the  depth  Vq 
which  yields  the  proper  binding  energy.   Bransden  has 
chosen  wave  functions  of  resonating  group  form  to  derive 
Integrodlf ferentlal  equations  for  the  radial  functions, 
and  has  solved  those  numerically  to  determine  the  scatter- 
ing phases.   He  has  found  a  potential  which,  assuming  near- 
Berber  exchange,  gives  the  right  phases  for  S  and  P  n-He 
elastic  scattering  and  n-D  scattering,  and  Is  consistent 
with  the  binding  energies  of  the  deuteron  and  the  alpha 
particle,  but  gives  too  small  binding  of  He'^  and  T,  and 
Is  of  too  long  range  to  give  the  best  fit  with  the 
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low-energy  two-body  scattering  data.   Burke  and 
Robertson  have  solved  the  deuteron  problem  for  a  gaussian 
potential  with  various  values  of  the  depth  V^^  and  the 

range  (x.   However,  Laskar  states  that  the  binding  energy 

7  R 
of  the  trlton  Is  8.5  mev,  while  It  actually  Is    8.483 

mev.   It  would  be  more  correct  to  say  that  this  is  a 
model  calculation.   The  fact  that  results  obtained 
with  the  polarized  gaussian  wave  function  agree  very  well 
with  experiment  must  be  viewed  as  a  fortuitous  coincidence. 

Instead  of  comparing  with  experimental  binding 
energies,  one  should  compare  results  with  various  theoreti- 
cal approaches.   The  above  model  of  the  trlton  was  solved 

9 
exactly  by  Baker  et  al . -^  for  various  values  of  the  parameters 

of  the  potential.   Kalos    has  solved  this  by  a  Monte  Carlo 

method.   Tang  et  al .    have  used  variational  principles, 

and  have  obtained  rough  agreement  with  the  abo\e methods . 

■:The  alpha  particle  with  no  Coulomb  forces  was  solved 

by  Kalos.    These  can  be  taken  as  a  guide  as  to  what  the 

true  energy  is,  to  see  how  accurately  the  polarized  wave 

function  gives  the  energy.   It  turns  out  to  lie  in  the 

range  given  by  the  above  calculations. 

Laskar  has  used  a  one-parameter  gaussian  trial  wave 

function  in  the  Rayleigh-Rltz  minimal  principle.   The 

parameters  of  the  potential  were  V^  =  -37-44  mev,  and 

[1  =  0.2669  fermi   .   He  has  obtained  an  energy  of  -6.744 

mev.   Comparing  with  the  experimental  value  of  -8.5  mev. 


we  see  that  the  simple  gausslan  is  a  very  poor  trial 
wave  function. 

3.   The  Polarized  Wave  Function. 

A  very  much  improved  wave  function  can  be  found 

by  considering  the  following  polarized  wave  function. 

The  polarized  function  can  be  pictured  as  follows.   We 

write  the  wave  function  (j)  =  >   cj) .  ,  where  the  functions 

i=l   ^ 
<\>.    are  functions  with  particle  i  in  an  asymmetrical 

position.   We  can  picture  the  situation  described  by  (j)^ 

as  follows.   We  consider  a  fixed  vector  x  joining 

particles  1  and  2,    and  a  vector  y  joining  particle  5  with 

the  center  of  mass  of  the  first  two  particles.   The  picture 

consists  of  y  being  perpendicular  to  x,  and  the  vectors 

having  average  values  for  magnitudes .   Although  the  actual 

case  consists  of  the  particles  arranging  themselves  so 

that  the  vectors  x  and  y  point  in  all  directions  and 

have  magnitudes  ranging  from  zero  to  infinity,  the  above 

picture  is  a  type  of  average  picture  convenient  for 

visualization  of  the  problem.   The  unpolarized  wave 

function  describes  a  case  where  the  three  particles 

lie  on  the  vertex  of  an  equilateral  triangle,  so  that 

2    H-      2 
X  =  -^  y   .A  polarized  wave  function  describes  the 

case  where  the  third  particle  is  further  or  closer  to  the 

center  of  mass  than  the  vertex  of  the  equilateral  triangle. 

The  polarized  wave  function  considered  here  is  a  linear 


10   - 


combination  of  the  two  types  of  polarizations. 

The  coordinates  are  defined  as  follows.   Let  r. 
be  the  position  vector  of  the  i-th  particle  and  R  the 

position  of  the  center  of  mass. 

1  ^" 

(1)  R  =  —  y        v.,    where  n  is  the  number  of  particles. 

^  i=l   ^ 
We  define  x  and  y  as  follows: 

(2)  X  =  Hi  "  II2 

(5)   1  =  £3  - 1  (i:i+Z2^- 

For  the  alpha  particle  we  have  in  addition  a 
coordinate  _z: 
(^)   z  =  r  1^  -  ^  (r^  +  r^  +  r^) 

We  define  the  square  radius  as  follows: 

(5)  r^  =1^  (r.  -  R)2  . 

^  i=l   ^ 
For  n  =  3  we  have 

(6)  r  =  -^  (x  +  ^  y  )  , 
and  for  n  =  4  we  have 

(7)  r^  =  i  (x^  +1  y2  +1  z2)  _ 

Equations  (6)  and  (7)  are  easily  proven  from  the 
definitions.   They  can  also  be  seen  geometrically  as 

follows.   For  an  equilateral  triangle  of  side  x,  height  y, 

4   2    2 
we  have  ^^  y  =  x  .   Let  the  distance  from  the  centroid  to 

2     2        '-'    1  4   2 
a  vertex  be  r.   We  have  r  =  x  /3  or  r^  =  ^(^  y  ) •   We 

take  the  average  of  these  expressions  to  get  equation  (6). 

To  derive  equation  (7)^  we  coinsider  a  tetrahedron  instead 

of  an  equilateral  triangle  and  proceed  in  a  similar  fashion 
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The  trial  wave  function  for  the  trlton  which 
Laskar  has  used  Is 

(8)  ct,  =  e-5r^^   , 

with  7  being  the  variational  parameter.   This  can  be 
written  as  follows.   Let 

1=1  ^ 

with 

-  |(x^ay2) 

(j).,  and  <|)p  are  defined  so  that  cj)  Is  symmetric  In  r.. 
For  a  =  4/3, we  have  Laskar 's  simple  gausslan.   For 
a  /  4/5,  we  have  a  polarized  wave  function.   The 

picture  described  by  ^^   is  as  follows.   ¥e  first  note 

2 
that  the  expectation  value  of  x   for  a  wavefunction 

7vx2 
e     is  1/A.   For  a  <  4/3,  therefore,  the  third 

particle  is  further  from  the  center  of  mass  than  the 

vertex  of  the  equilateral  triangle,  while  for  a  >  4/3, 

the  third  particle  is  closer.   (j)^  and  cji^  describe  the 

cases  where  particles  1  and  2  respectively  are  displaced 

from  the  vertex. 

We  proceed  as  follows.   We  first  fix  a  at  4/3,  then 
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vary  <H>  as  a  function  of  7  until  we  find  a  minimum.   We 
fix  7  at  this  value,  and  write 

-  ^(x  +ay  ) 

and  vary  f  and  a.   f  has  the  meaning  of  a  scale  factor, 
since  If ;«  push  one  particle  out,  we  expect  to  reduce  the 
size  of  the  system,  and  vice  versa.   There  are  two 
geometrically  different  types  of  polarizations: 

4 
a  <  T^  ,   f  <  1     (extended) 

4 
a  >  ^  ,   f  >  1    (compressed) 

Either  case  gives  roughly  the  same  Improvement  to  the 
energy,  with  the  extended  case  slightly  better.   A 
linear  combination  of  these  leads  to  a  resonance  effect, 
giving  a  large  Improvement  to  the  energy.   This  wave 
function  Involves  five  parameters.   It  Is  written  as 

2 


(j)^  =  >     C^  ^-^. 

(9) 


-  2f^^^  +a^y") 

c^  Is  determined  by  normalization  to  be  equal  to  one 
For  Laskar's  potential,  we  obtain 
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^^2y/2  ^  2^55  f  E  =  -  8.496  mev 

f-L  =  0.4131       J  ^i  =   0.3374 
c^  =  0.146?         r^  =   5-2471       I  ^2  =  2.096 


ci  =  1. 


7  was  set  equal  to  .2357?  which  is  Laskar's  value. 

7  8 
The  experimental  energy  is  E  =  -8.483  mev,  ' 

and  the  root  mean  square  radius  as  seen  by  electron 

1 2 
scattering    is  1.68  +  0.I7  f. 

The  unpolarized  wavefunctlon  gives  the  following 

E  =  -  6.67  m.ev.         <r^>^^^   -  1.48  f. 

We  see  that  the  polarization  causes  the  radius  of  the 
triton  to  increase. 

Instead  of  restricting  a.  and  f.  so  that  the  wave 
function  can  be  interpreted  as  a  resonating  sum  of 
polarized  functions,  the  parameters  can  be  allowed  to 
vary  over  various  arbitrary  ranges.   Many  sets  of  values 
of  the  parameters  are  found  which  yield  minima  f or  <  Hj>. 
This  is  to  be  expected,  as  a  many  body  function  can  be 
approximated  in  many  ways.   However,  none  of  these  minima 
are  as  low  as  the  value  found  by  the  physically  meaningful 
restriction . 

If  one  adds  a  gaussian  term  to  the  polarized  wave 
function  so  as  to  get  a  function  involving  7  parameters, 
one  finds  a  negligible  improvement. 

For  the  alpha  particle,  we  proceed  in  an  analogous 
fashion.   We  replace  r^  by  ■i(x^+  ay^+  bz^)  .   The  picture 
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for  the  unpolarized  case  Is  that  of  four  particles 
situated  at  the  corners  of  a  tetrahedron.   There  are 
three  possible  types  of  polarization.   One  consists  of 
two  particles  closer  to  the  center  of  mass  than  the 
vertices,  and  the  second  case  is  with  two  particles 
further  away.   The  third  case  is  with  one  particle 
closer  and  the  other  particle  further.   That  is, 
either  a  >  4/3,  b  >  3/2,  f  >  1,  or  a  <  4/5,  b  <  3/2, 
f  <  1,  or,  for  the  third  case,  a  <  4/3,  b  >  3/2.   A     ^ 
fourth  case,  a  >  4/3,  b  <  3/2,  does  not  occur  as  the 
particles  are  identical,  this  case  being  identical  to 
the  third  case.   A  trial  function  consisting  of  a 
linear  combination  of  the  first  two  types  of  polariza- 
tions yields  excellent  results .   Addition  of  the  third 
case  gives  an  order  of  magnitude  less  Improvement,  and 
so  to  a  good  approximation  the  wave  function  of  the  alpha 
particle  can  be  considered  as  the  sum  of  two  polarized 
functions . 

For  a  simple  gaussian,  i.e.,  the  unpolarized  wave 
function,  the  minimum  in  the  energy  can  be  found 

analytically.   This  serves  as  a  check  on  the  calculations, 

g 
For  the  potential  of  Baker  et  al.-^,  Vq  =  -5I.5O  mev.  and 

[1  =  1/(1. 6f)  ,  the  minimum  is  at  E  =  -6.803  mev. 

Tang  et  al .    get  a  value  of  E  =  -6.9I  mev.   This  is 

?  ?  2 

because  they  set  -fi  /m  =  41.47  mev/f  .   I  have  set  -h  /m 

2  2 

=  4l,59  mev/f  .   In  computing  -R  /m,  the  value  of  -fi  was 
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13 

taken  from  a  recent  reevaluation  ^  of  fi.   The  value  of 

the  mass  was  computed  from  the  miass  defects  given  in 

14 
Lelghton   .   Actually,  there  is  no  need  to  be  this 

precise  about  the  mass,  for  relatlvistlc  effects  are  not 

considered  here.   For  the  triton,  relatlvistlc  effects 

are  about  .7  mev.,  as  shown  in  the  appendix.   For 

Laskar's  potential,  the  minimum  is  at  E  =  -6.675  mev. 

Laskar  uses  Tang's  value  of  -h  /m,  and  gets  an  energy 

of  E  =  -6.744  mev. 

The  analytical  expression  for  the  minimum  of  <H^ 
is  a  fifth  degree  equation,  which  has  two  real  positive 
roots.   The  second  root,  however,  gives  a  positive  value 
for  the  energy,  indicating  that  there  are  no  excited 
bound  states.   The  energy  here  is  E  =  +.8  mev. 

The  following  results  have  been  obtained  for  the 
model  trlton.   Baker  et  al .  have  obtained  E  =  -9-47  mev. 
Kalos   has  used  the  same  potential  and  has  obtained 
E  =  -9.42  mev.   Tang  et  al .  have  obtained  E  =  -9-74 
mev.   The  polarized  wave  function  yields  an  energy  of 
E  =  -9.52  mev.  with  fi^m  =  41.59  mev/f^,  and  E  -  -9.62 
with  fi^/m  =  41.47  mev/f^.   It  is  interesting  to  note 
that  Tang  et  al .  have  obtained  almost  this  same  value 
of  the  energy,  E  =  -9.63  +  0.04  m.ev,  with  a  gaussian 

type  of  wave  function: 

2  2 

3     -arf.      -pr.. 

(j)  =  -pi   [e   ^J  +  c  e   ^'^l 

i<j=l 
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which  involves  three  parameters  instead  of  five.   It  is 
quite  possible  that  this  function  is  very  similar  to  the 
polarized  gaussian  function,  for  a  function  of  many 
variables  can  usually  be  written  in  many  ways.   However, 
Tang's  form  does  not  permit  an  easy  generalization,  as 
it  has  no  simple  physical  interpretation.   It  appears 
that  the  polarized  gaussian  wave  function  is  the  best 
simple  function  for  the  energy. 

Another  check  on  the  calculations  is  found  in  the 
virial  theorem.   The  following  theorem  is  proven  in  the 
appendix:   If  a  variational  function  permits  a  scale 
transformation,  the  value  of  the  variational  parameter 
which  gives  a  minimum  of  the  energy  also  satisfies  the 
virial  theorem  exactly.   The  wave  functions  which 
minimized  the  energy  actually  did  satisfy  the  virial 
theorem,  to  machine  accuracy.   However,  first  order 
changes  about  the  minimum  produced  first  order  (or 

larger)  changes  in  the  violation  of  the  virial  theorem. 

3       4 
For  He  and  He  we  have  to  take  into  account  the 

Coulomb  potential.   Since  this  is  small,  an  approximate 

expression  for  the  Coulomb  potential  was  substituted. 

For  a  simple  gaussian  wave  function,  the  Coulomb  energy 

can  be  calculated  exactly.   For  small  changes  in  the 

wave  function,  the  Coulomb  energy  will  change  by  a 

small  amount,  which  is  approximately  inversely  proportional 

to  the  root  mean  square  radius.   Consequently,  the  Coulomb 

energy  was  set  equal  to  a  constant  divided  by  the  root 
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mean  square  radius,  and  the  constant  was  determined  from 
the  unpolarized  vjave  function. 

The  minima  were  found  by  a  random  function  minimization 
routine,  called  L0C^,  developed  by  Jay  Leavitt.   One  gives  a 
fixed  region  in  parameter  space,  and  the  routine  searches 
randomly  over  this  space,  converging  about  a  minimum. 
It  contains  features  which  help  it  to  overlook  local 
minima.   Another  routine  used  was  SC00P,  developed  by 
Mel  Sobol.   This  routine  varies  the  parameters  singly, 
then  simultaneously,  in  a  manner  avoidng  the  pitfalls  of 
the  ordinary  gradient  methods.   It  is  good  for  a  function 
of  three  or  less  parameters,  but  for  more  parameters  L(2('C^ 
is  necessary.   SC0i2^P  has  been  submitted  to  SHARE,  whereas 
L^C0  has  not . 

The  success  with  the  polarized  gaussian  potential 
suggests  that  a  polarized  wave  function  be  used  in  nuclear 
physics  problems,  instead  of  a  wave  function  which  is  a 
product  of  two-body  functions,  or  a  function  of  a  quantity 
invariant  under  particle  exchange.   A  polarized  function 
would  be  about  as  easy  to  handle  as  an  unpolarized  function, 
and  will  probably  give  much  better  results.   One  introduces 
the  polarization  as  follows.   First  one  replaces  quantities 
invariant  under  exchange  of  particles  by  polarized  quantities 
For  a  three  particle  system,  the  polarized  quantity  merely 
distinguishes  the  third  particle.   For  example,  instead  of 
a  sum  of  the  particle  distances,  we  would  write  this  sum 
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with  the  distance  of  the  third  particle  multiplied  by  a 
parameter.   If  the  three  particles  are  Identical,  we 
replace  our  original  function  Of  the  Invariant  quantity 
by  a  sum  of  functions  with  the  1-th  particle  multiplied 
by  the  parameter.   This  gives  us  a  polarized  function. 
We  then  Introduce  resonance  by  replacing  our  polarized 
function  by  a  linear  combination  of  polarized  functions, 
with  the  parameter  greater  and  smaller  than  one.   For  a 
four  particle  system,  it  probably  would  be  good  enough 
to  multiply  the  third  and  fourth  particles  by  the  same 
parameter,  and  then  take  a  linear  combination  of  two  types 
of  polarizations.   This  is  because  in  the  case  of  a 
gaussian  potential  and  a  polarized  gausslan  wave  function, 
the  third  type  of  polarization  gives  us  a  small  improve- 
ment . 

For  He   it  may  be  possible  to  give  different  polariza- 
tions to  the  neutron  and  protons,  permitting  a  more  direct 
calculation  of  the  Coulomb  effects  and  of  charge  independence 
of  nuclear  forces. 

Efimov   has  used  a  type  of  polarized  gausslan  wave 
function  for  a  calculation  of  the  scattering  lengths  for 
elastic  scattering  of  neutrons  by  deuterons.   Massey   has 
discussed  this  work  and  similar  work.   Eflmor's  wave 
function  of  the  system  neutron  plus  deuteron  is 
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i^{j,r)    =  -  e'^y   -  y^   Cp^  exp  [-y^[ix^   -   nuy^] 


-  JZ  ^m  exp  [-m^L(x^+y^)/2]  , 
m 

with  x  the  relative  coordinate  of  neutron  and  proton 
within  the  deuteron,  and  y  that  of  the  Incident  neutron 
relative  to  the  center  of  mass  of  the  neutron,  and  with 

c   ,  A  adjustable  parameters,   u-  determines  the  range 

2 
of  the  potential  e~^   ,  and  7   are  quantities  occurring 

P 

In  an  approximation  to  the  deuteron  wave  function  ^{x) 

^  2 

Mx)  ^^I  b^  e  P 
p=l  P 

A  greater  effect  was  obtained  with  the  wave  function  of 

Eflmov   than  with  the  other  types  of  polarized  wave 

functions  discussed  by  Massey.   Eflmov 's  wave  function 

resembles  the  polarized  wave  function  of  this  paper,  but 

does  not  Incorporate  the  ideas  as  much.   It  is  probable 

that  a  calculation  done  with  this  type  of  polarized 

wave  function  would  yield  much  better  results. 

4 .   The  Three-Body  Correction  Potential. 

Let  tjj   be  the  true  wave  function,  which  satisfies 

(E-H)^  =  0  , 
and  let  cf)  be  the  polarized  wave  function  which  minimizes 
the  energy.   ¥e  would  like  to  know  more  about  (f)  and  how 
it  differs  from  the  true  wave  function  ^.   We  can  start 
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In  that  direction  by  the  introduction  of  a  three-body 
correction  potential  W  which  is  defined  so  that  (E-H)<t)  =  Wcj), 
i.e.,  so  that  our  trial  wave  function  (j)  is  an  exact  eigen- 
function  of  the  hamiltonian  formed  by  adding  W  to  the 
original  hamiltonian.   W  can  be  expressed  in  various 
ways.   One  way  is  by  using  projection  operators.   Let 
P  be  the  projection  operator  |  (j)>  <C  (j)  |  .   We  have 

(10)  -  W  =  PqA(I-Pq)  +  (I-Pq)HPq  . 

'We  easily   see  that  -W({)  =  (H-E)(t)  and  -([)W  =  (|)(H-E)  . 

Another  way  to  express  the  three  body  correction 
potential  W  is  in  configuration  space.   We  define  W  as 
the  operator 

(11)  W(r^,r2,r:^)  =  -^^^   • 

W  appears  as  a  sum  of  gaussians  divided  by  a  sum  of 
gaussians.   It  has  the  form  of  a  function  which  becomes 
large  near  Snd  very  far  from  the  origin,  with  two  minima 
around  the  mean  square  radius,  due  to  the  polarizations 
in  the  wave  function.   It  does  not  start  becoming  large 
again  until  the  wave  function  becomes  very  small. 

We  could  find  a  potential  W„  which  behaves  like  a 

3. 

physical  potential,  i.e.,  goes  to  zero  at  infinity,  yet 
approximates  W  in  the  regions  where  the  wave  function  is 
not  small.   To  begin,  we  note  that  (W^  is  almost  zero, 
and  so  W  does  not  contribute  to  the  energy.   This  follows 
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from  the  fact  that  <(H  +  W  -  E>  =  0  and  <(H  -  E/>  Is  almost 

zero,  as  the  polarized  wave  function  Is  a  good  wave  function 

for  the  energy.   We  therefore  require  that  < W  >%  0.   We 

a 

cannot  write  the  wave  function  as  a  simple  gausslan  times 
an  expansion  In  x  and  y,  and  use  this  to  find  W„ ,  as  the 

a, 

polarization  of  the  wave  function  Is  too  large  to  permit 
an  expansion.   We  approach  It  as  follows.   We  set  W„ 

3. 

equal  to 

2 
-Zl^        2     4     4 
(12)  e      (q+Zgr  +z  r  +  Z|^s  )  , 

2 
where  q  Is  the  value  of  W  at  x  =  y  =  0,  and  r   is  the 

4 
square  radius  defined  In  equation  (1),  and  s   Is  the 

polarization  radius,  defined  as  follows.   Let 


(13)  2^3  =  Zi  -  1:2  ' 

(14)  13  =  1:5  -  ^(£1+1:2^  ' 


and 


similarly  define  x-i  >  xp'  Zi  ^  ^^'^   Zp-   ^^  have 


(15)  s^  =^  (xi  •  Zi)^  • 

The  z.  are  parameters.   The  factor  multiplying  the  exponent 
Is  the  simplest  fourth  order  polynomial  syrrimetric  In  the 
particle  distances.   To  find  the  parameters  z.,    the  follow- 
ing expression  was  minimized 

^[(«a  1,3,1c  -Wl,J,l.'*i,J,l^''+l<"a>l  • 

2    2 

W,  W^,  and  (})  are  functions  of  the  three  variables  x  ,  y 


22 


and  X  '  Z  •   These  functions  were  evaluated  for  many 
values  of  the  variables.   The  subscript  (i,j,k)  indicates 
that  the  function  is  to  be  evaluated  for  a  given  set  of 
parameters.   We  get  W„  =  W  to  lO/o  accuracy,  and  have 

a. 

(^^y  =    2  X   lO"-'  mev.   W^^  begins  to  differ  from  W  in 

regions  where  (})   =  .01  ^{0)    ,    where  ^{0)    is  the  value 

of  the  wave  function  at  the  origin.   W^  can  be  used  in 

a 

expressions  involving  expectation  values  of  quantities 

which  are  small  at  large  distances  from  the  origin. 

This  is  because  (j)   is  extremely  small  in  the  regions 

where  ¥^  is  very  different  from  W. 
a 

Using  Laskar's  potential,  the  parameters  obtained 

for  W  are 
a 

^-^   =  .5756,   q  =  17.153,   2^  =  4.154,   z^  =  -2.541, 

zj^  =  11.053. 

One  can  try  many  different  expressions  for  W_ 

a. 

involving  more  parameters .   A  negligible  improvement 
is  obtained. 

This  approximate  expression,  W  ,  has  the  following 
utility.   If  the  wave  function  of  the  triton  were  actually 
a  polarized  gaussian  of  the  form  found  here,  then  there 
would  exist  a  three  body  potential,  which  can  be 
approximated  by  W^ .   One  can  proceed  to  solve  the  actual 

a. 

problem,  the  three  body  nucleus,  with  a  two  body  gaussian 
potential,  as  follows.   Write 
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where  Hp^  =  H  +  W  .   We  expand  in  a  perturbation  series, 

knowing  that  H^cj)  %  E(}) .   We  know  that  W  is  small  --  its 

expectation  with  respect  to  (j)  is  zero.   W^  is  largest  at 

a 

the  origin,  when  it  is  about  17  mev .   H(t)/(j)  is  about-26  mev 
here.   Although  W  is  very  large  at  the  origin,  it  falls 
off  very  fast,  about  as  fast  as  (j). 

W„  does  not  have  the  proper  asymptotic  behavior. 

a. 

Consequently,  it  will  not  give  correct  results  when 

values  of  W^  are  needed  for  large  r.   For  example,  for 
a 

Laskar's  potential  we  have 

(^^}  =   2.84?  mev^ 
and 

(w^>  =  8.743  mev^ 

The  ratio,  (v!^}/(^^},    is  5 -06. 

The  three  body  correction  potential  can  be  used  to 
study  the  second  order  change  in  energy.   We  write 

-  AE  =  ZI  Y^ 
n>0  ^n  ^0 

2 


1_  lw 


-  E,  -E^  ^-rr     '  nO 

1   0  n>0 


Now 


y~  IW^^I^  =YI1    (Olw|n)(niw|0) 
n>0   ^^^     n>0  ^ 


=  XZ  (0|w|n)  (n|w|0),   as  (•0|w|0)=0 
n>0  - 

=  (o|w^|o)  . 
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We  therefore  get 

(16)  _AE<^^<w2>. 

For  Laskar^  potential  we  get,  setting  E^  =  -  .4  mev, 

-  AE  <  -^3^  =  1.08  mev. 

For  the  potential  of  Baker  et  al .  we  get,  where 
W  /  =  16.9^  niev  ,  and  E^  -   -A   mev,  as  discussed  below 

-  AE  <  ^4^  =1.86  mev. 

The  method  of  Temple  (see  Schmid  et  al .    )  consists 
of  varying  parameters  to  maximize  E+AE.   The  result,  E+AE, 
Is  a  lower  bound  on  the  energy.   The  formula  for  the  lower 
bound  is 

r2\     y,,v2 


(17)  E>E^=<H>- 


<H^>  -  <H>' 
~  E^  -<H> 

subject  to  the  condition  that  <'H^  is  less  than  E-,  . 
E^  Is  the  eigenvalue  of  the  first  excited  state  with  the 
same  symmetry  property  as  the  ground  state.   This  is 
taken  as  the  energy  of  the  deuteron.   We  assume  that 
there  are  no  excited  bound  states,  as  indicated  by 

calculations  with  a  gaussian  wave  function.   The  value 

11a 
of  E,  was  evaluated  by  Tang  et  al.    ,  and  is  equal  to 

-0.40  mev  for  their  potential.   In  this  expression, 

<H  >  is  given  simply  by  <^(T+V)  >  ,  as  spin  forces  are 

not  considered.   The  result  for  E^  is  E^  =  -10. 69  mev. 
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5.   An  Upper  Bound  to  the  Error  In  the  Mean  Square  Radius. 

A .    General  Calculations. 

The  polarized  gausslan  wave  function  has  been  found 
to  be  an  excellent  approximation  to  the  true  wave  function. 
It  is  of  interest  to  see  how  good  this  wave  function  is 
for  physical  quantities  other  than  the  energy.   At  first 
sight,  it  appears  that  this  wave  function  is  not  good  at 
large  distances.   To  begin,  the  true  wave  function  behaves 
asymptotically  like  an  exponential,  rather  than  a  gausslan. 
Furthermore,  the  three  body  correction  potential  W  becomes 

very  large  at  large  distances,  diverging  at  infinity  like 

p 
r  ,  suggesting  that  the  polarized  wave  function  is  not 

correct  at  large  distances.   However,  the  wave  function 

becomes  very  small  before  W  starts  becoming  large.   To 

study  the  question  more  fully,  an  upper  bound  on  the 

mean  square  radius  is  calculated  in  this  section. 

The  result  turns  out  to  be  very  large.  This  means 

that  it  is  possible  that  the  true  wave  function  gives  a 

result  for  the  mean  square  radius  very  different  from 

that  calculated  with  the  polarized  wave  function. 

However,  it  does  not  mean  that  it  will  necessarily  be 

very  far  off.   To  get  a  better  view  of  what  actually 

happens,  a  variational  principle  for  the  mean  square 

radius  is  studied  in  Section  6,  and  there  the  result  is 

very  good.   However,  the  following  upper  bound  calculation 
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Is  of  Interest  in  Itself,  as  the  method  can  be  applied 
to  other  areas  of  physics. 

To  calculate  an  upper  bound,  one  may  try  to  apply 

17 
the  results  of  Redel.  '  Redel  has  derived  an  upper  bound 

for  the  error  in  the  density.   His  upper  bound,  however, 

contains  a  constant  term.   If  one  writes  the  error  in 

the  mean  square  radius  as  bounded  by  the  integral  of  the 

mean  square  radius  times  the  bound  on  the  density,  one 

will  get  an  upper  bound  of  infinity  due  to  the  Integral 

of  the  mean  square  radius  times  a  constant.   Hence 

Redei's  method  will  work  only  for  finding  an  upper  bound 

on  the  error  in  some  quantity  which  goes  to  zero  at 

infinity  fast  enough.   Consequently,  another  approach 

for  finding  a  bound  on  the  error  in  the  mean  square 

radius  has  been  derived,  as  follows. 

Let  f   be  the  correct  wave  function,  i.e.,  (E-H)?//  =  0, 
and  let  (f)  be  the  polarized  gaussian  wave  function  found  by 
minimizing <  H >.   We  want  to  get  a  bound  on  the  difference 
between  (^,r  ^)  and  {^,r   <\))  .      We  proceed  as  follows. 

In  order  to  evaluate  (^,r  ?//),  we  express  the  trial 
wave  function  ^   in  terms  of  ip   and  the  Green's  function  G, 
where  G  is  the  principal  part  of  (E-H)~~,  as 

(j)  =  ^  +  GWcj)  . 

We  normalize  (^,^)  =  1.   Let  P  be  the  projection  operator 
to  the  state  tp ,    i.e.,  P  =  li'^^Cj^l,    and  let  Q  =  1-P .   Q,  is 


27  - 


(19a) 


the  projection  operator  to  the  complement  space. 
We  have 

((j),^)  =  (^,^)  +  (^,GW(t))  . 

G  can  be  defined  as  G  =  Q(E-H)~  Q.   Consequently,  the 
second  term  In  the  above  equation  vanishes,  and  we 
have  ((f),^)  =  1. 

Summarizing,  we  write 

(l8a)  (E-H)^  =  0 

(l8b)  (V/,^)   =  1 

(E-H)(j)   =   Wcj) 
^   =   Tp   +  GWcf) 

(19b)  ((f),^)    =    1 

(20)  P   =     |V/><^| 

(21)  Q  =   1    -   P 
(22a)  G  =  QG^Q 
(22b)  G^    =    (E-H)"-"-    . 

We  can  get  an  upper  bound  on  ((|),(1))  as  follows. 
Define  a  normalized  function  f  orthogonal  to  f,    and  a 
quantity  r| ,  giving  a  measure  of  the  error  In  (j): 

(23a)  ^  ^,^  =     yT^  ^  +  Tif  , 

(9.9) 

(23b)  (^,f)  =  0  ,     (f,f)  =  1  . 
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If  we  let  E-j^  be  the  energy  of  the  first  exlcted  state, 

as  above,  we  have  the  Inequality,  assuming  the  ground  state 

of  H  is  non-degenerate: 


Hi 


172  -   '/^  E^  +  ^%^ 


Using  this,  we  have  for  an  inequality  for  t^ 

((t),H(t)) 


>  (l-Ti^)E  +  ri^E^ 


Let  Ey  =  ((t),H(j))/((j),c{)) .  E„  is  an  upper  bound  on  the  energy, 


(24) 


%  1  E  +  ri^(E^-E) 


2   ^0"^ 


'1 

To  find  E,  we  can  look  for  a  lower  bound  on  the  energy,  E^. , 
as  discussed  above  in  Section  4.   t]  has  the  following 
significance.   We  write  cj)  as 


Comparing  with  (19)^  we  see  that  'y{<^,^){l-r\    )  =  1,  as 

we  have  in  both  cases  written  cj)  as  the  sum  of  two  orthogonal 

functions.   We  then  have 

(25)  i^A)  =  — ^  • 

l-Tl'^ 

2 
We  can  also  find  an  upper  bound  on  t)   as  follows. 

We  define  expectation  values  with  respect  to  (j),  i.e., 
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^9/>  =  ((j),0(())/((j),(t))  for  any  operator  0. 

(^,(t))  =  1  +  <WG^W>(({),(1))  . 

Solving  for  ((jj^^))  and  comparing  with  equation  (25),  we 
see  that 

¥e  find  an  upper  bound  as  follows-   Let  /t  be  any  function 
satisfying  the  conditions  ip   satisfies.   Define  a  number  G^^ 


so  that 
(26a) 
We  have 


(^,G^) 


(;^,;^) 


^%  • 


We  set 
(26b) 

We  have,  therefore. 


%^ 

1 

^-^1 

• 

^0   = 

1 

• 

^1-^ 

(27) 


Ti^   <  g2<w2> 


,2  /,,2' 


It    turns    out    that    (Ejj-E)/(E^-E)    <  Gq  <W  ^   . 

2  '  '  y        2     \ 

Consequently,  -r]      can  be  used  as  an  upper  bound  on  <_WG  W^, 

2 
where  we  set  i^   equal  to  the  upper  bound: 

2  _  %"^ 
We  can  try  values  of  E  lying  In  the  range  E^  <  E  <  E^  . 
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Define 


(28)  AZ  =  (^,Z^)  -  <Z> 


for  any  quantity  Z.   AZ  Is  the  error  In  the  expectation 

value  of  Z.   We  are  Interested  In  finding  an  upper  bound 

I   2i 
for  |Ar  |.   This  can  be  expressed  as  follows. 

Ar^  =  I(^,r2^)  -  ((t),r2(t,)  +  ((t),r2(f))  [1  "  j^l  I 
Ar^  _^  \{ij,r^ij)    -    (({),r^(l))  I  +  r]^{^,r^^)     . 


Judging  from  atomic  physics  problems,  E  Is  probably 
a  lot  closer  to  K,  than  E- .   Schmld  et  al .    show  that 
the  same  Is  true  for  the  deuteron  bound  by  a  gausslan 


potential . 

We  want  an  upper  bound  on  |  (ij!/,r'^?//) -(({), r'^cf)) 

(29)  \if,r^ij)-{^,r^<\>)\    =  |(^-i>,r^(^-(t)))  +  2  Re{f-^,r^^)  \ 

1 

(^-^,r^(^-(t)))  +  2  Re((^-<t)),r^(^-(t)))^ 

1- 
X  (<{),r"(j)) 

In  the  above  expression,  either  (■<p-^,v'^^)    >   0,    or  else 
we  write  ^-fp   throughout. 

(30)  I  (^,r^^)-((t),r2(}))  I  <  (^-(t),r^(^-ct))) 

+  2  Re(^-4),r^(^-^))^/'^<(r^^/2 
1 


V^ 
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(31)  (^-(l),r^(^-<j)))  =  ((l),(l))<WCTr^GW>  . 

From  equation  (26)  we  know  an  upper  bound  on  the  diagonal 
elements  of  G.   This  upper  bound  cannot,  however,  be 
substituted  in  equation  (31)>  as  there  the  off  diagonal 
matrix  elements  occur.   It  is  necessary  to  commute  r  with 
G  to  get  an  expression  like  equation  (26).   This  has  been 

done  in  the  appendix.   The  result  is,  where 

I  X       1 

r.  .  =  r .  .  -  -7^  r,  .  -  -?5  r .  .,  i,k,i  cyclic, 
ij     ij    2   kj    2  ij'   '  ' 

(32)  <¥Gr^GW>  =  <(wrG^rW>  +  ^  ^  Im  ^^  <WGp .  .G^r'  .¥> 


+  1    (|)'lZ<WGPijG2p:.GW> 


We  apply  Schwartz's  inequality.   Note  that  equation  (a6i) 
from  the  appendix  is  {)C,v'^.y.)    =^{'y.  ,v^X^  >    where  X^   ' 
This  expression  can  be  evaluated  exactly.   Schwartz's 


inequality  applied  to  the  second  term  of  equation  (32)  is 
(33) 


4  ;R 
9  m 


Im  "^^  <^¥Gp .  .Gr .'  .¥>  | 
<  4|  ||Gp.jGW(|)||    |lGr^jW(t)||/(()),<j)) 

=   4G2[(^)<¥r2w>]l/2  ^    ^ 
where  Kq    is   an  upper  bound   on   the   kinetic   energy,    defined   as 


32      - 


(GW(}),KGW(j))  <   KQ(GW^,GW(t))  . 

Kp.   will  be   evaluated  below.   k  =  1i  /2m  . 

In  applying  Schwartz's  Inequality  to  the  last  term  of 
equation  (52),  we  note  that  equation  (A62)  is 


We  also  have 


and  so 


(  X  ,P^^X  )  =  9(/t^P  />^) 


l|Gp'  .GW(1)1|        

— Ifi^ —  <  Go  v^^  i  ^ 


We  get  for  our  term: 
2  f-fis^ 


(3M  If  (h)     IZ  <wgp..gS..gw> 


2         1       ,o   T^    \   ^^   r>2      2 


4  2      '^ 

=   27   ^0  ^      Gq 


We   finally  get 


(35)     <(WGr^GW>  <  GQ<Wr^W>  +  4G^[-^<Wr^W>]-'-/^  ^   +  2?  ^o^^^O 


Setting  the  right  hand  side  equal  to  GRG,  we  have  as  our 
final  result  for  the  error  in  the  mean  square  radius: 

(36)  |(^,r2^)-<r2>|  <~^   [GRG+2[GRG<r2>]^/^  +  ^\r^>]     . 

1-Ti^ 
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The  upper  bound  turns  out  to  be  several  times  as 
large  as  the  mean  square  radius.   This  is  due  to  the 
large  value  of  -q,  or  ^H  />.   Although  ^E-H>  is  very  sm.all, 

■P  \  'P 

-<^(E-H)  y  is  not  small,  being  instead  of  the  order  of  E  . 
However,  in  problems,  such  as  atom.ic  physics,  where  t]  is 
small,  this  method  should  work,  and  give  a  good  upper 
bound  for  the  mean  square  radius. 

One  can  try  to  vary  the  parameters  of  the  polarized 
wave  function  so  as  to  minimize  the  upper  bound  on  the 
mean  square  radius.   The  improvement  is  very  small.   The 
upper  bound  is  much  larger  than  the  value  of  the  mean 
square  radius. 

Although  the  possible  error  on  the  mean  square 
radius  is  quite  large,  the  probable  error,  as  seen  by 
the  variational  principle  in  the  next  section,  is  very 
small . 

B.    An  Upper  Bound  for  an  Expectation  of  the 
Kinetic  Energy  • 

In  order  to  complete  the  above  calculation  on  an 
upper  bound  to  the  error  in  the  mean  square  radius,  it 
is  necessary  to  know  an  upper  bound  to  |lpG¥(l)|l  .   This 
can  be  expressed  in  terms  of  the  function  f  defined  in 
equations  (23) .   We  have 

(37)  ri(l-ri2)"^/^  f  =  G  W  (1)  . 
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The  procedure  Is  as  follows.   A  Coulomb  potential  V^ 

is  introduced  for  mathematical  reasons.   Generalizing 

1  P, 
the  work  of  Kinoshita   ,  we  can  find  an  upper  bound  for 

the  expectation  of  the  kinetic  energy  (f,Kf)  in  terms  of 

(f.v'^f).   Relating  (f,V^f)  with  (f,Vf),  where  V  is  our 

gaussian  potential,  and  expressing  (f,Hf)  in  terms  of  r|, 

we  have  our  upper  bound. 

Let  our  Coulomb  potential  V*^  be  defined  as 


where 


and 


0  ^—  '  ij 


c2  ^ 
^0     m 


-ij   -i   -j 

Let  X  =  r^p,  w  =  -2(Zi+I!2^  ^'^^  '^  ^    '— '  '  '^'^'~    '^^^'^''^^ 
element  dr  is  dx  dw.   Let  5  =  h/'dx. 

We  start  with  the  expectation  value  of  our  Coulomb 
potential,  then  integrate  by  parts  and  use  Schwartz's 
inequality.   Let  dQ  be  the  solid  angle  elem.ent  in  x 
space.   We  have 

(58)  (f,  I  f)  =ydw  dO^  dx  X  f2 

=  -  /  dw  do  dx  ^  S   f^  , 
integrating  by  parts.   We  rewrite  this  as 
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(f,  -  f)  =  -  /  dw  dO  dx  x^  f  S  f 
^  '  X        J   —    X  X 


(f.  ^   ^^  1   (/^^  ^^)^^^  [y^dT(S^f)2]l/2 


by  Schwartz's  inequality.   Since  f  is  normalized,  we  have 

(f,  ^  f)  <  [TdTfa  f)2]V2  _ 


Since  dr  and  f  are  symmetrical  in  the  particle  indices, 
we  have 

(59)      (f,V^f)  <  3  V^[  y^dT(S^f)2]l/2  . 

We  define  the  vectors  x. ,   ^.    like  in  Appendix  6,    where 

for  example,  x^  =  L'   K-a  =   Z'   '^^®  index  i  tells  us  which 

particle  is  to  be  distinguished.  Our  wave  functions  can 

be  considered  as  functions  of  the  variables  {]^^,J^,^)  • 


We  have 


Sf  ^  y—  Sf   .  ^-i 


Using  the  following  equations 


we  get 


1     1 

Z3  =  I  Z  +  R   ^ 


Bf   _  1  Sf    1  Sf 
"SxZ   2  "Sr^  "  2  "Sr^ 
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by  considering  f   a  function  of  (x^^Y^^R) • 
Considering  f  a  function  of  (x-,  ,  y.,  ,R )  we  get 


Bf  ^  1  Sf  _  1  Sf 

dx.    2  dr^   2  "StZ 
—1      — d      —5 


Likewise  we  get 


Sf  ^  1  Sf  _  1  Sf 

"Sxl   2  "StT   2  "SFT  • 
—d  — i       — ;) 

Summing  the  squares,  we  get 

.r —  /-Sf  x2  _  1  ^; —  ,Sf  s2    1  ^; —  Sf     Sf 
A-  ^^^^   -  2  ^^  ^^F7^   ~  2  <— ^  ^F7  '  ^FT  • 

1-1  1   .  -1         KJ   -1     -J 

We  use 

,Sf  ^2   .; —  ,Sf  n2  ^  „  ^; —  Bf    Sf    „ 

as  f  does  not  depend  upon  the  center  of  mass,  and  get 
^—  ,af  v2   l,Sf  2  _  3  5—  /Sf  n2 


—1       —        1—1 


or 


1     —1         1     —1 

The  function  f  does  not  depend  upon  the  six  components 
of  the  vectors  x  and  y,  but  only  on  three.   Let  these 
three  be  \x.\  ,     |y  |  ,  and  an  angle  giving  the  orientation 
of  the  vector  y.   We  then  have 

1      1 
This  Is  shown  as  follows. 
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dx.        "Sx  X  ' 

as  f  does  not   depend  upon  any  angles  associated  with 
the  vector  x-   Squaring  and  summing  gives  equation  (4l). 
Equation  (39)  becomes 

(h2)  (f,V^f)  <|v^  f  ^/^^  ^^f^^^^ 

or 

(^5)  (f,V^f)  <  ^-   (f,Kf)^/^  . 

¥e  want  to  fin<3  a  V,  such  that 

2 

-ur . 

(44)    -  V^(f,Vf)  <  (f,V^f)  ,    where  V  =  -V^Se   ^^ . 

This  is  satisfied  if  for  all  r 

2   V^ 

v^Vq  e     1    r     ' 

Asymptotically,  for  small  and  large  r,  this  is  satisfied 
for  any  V^  .   Consequently,  if  at  some  point  the  left  hand 
side  is  equal  to  the  right  and  the  derivatives  are  equal, 
the  inequality  will  be  satisfied  everywhere. 
Let     ^2  ^   ^'o/^l"^0  '   ^^  require 

2  ^  g-M,r 


2 

2^Lre  ^ 


We  get 
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and  so 

Solving  for  V-,  ,  we  have 


or 


(^5) 


-l- 

2M,e  r? 

,2 
1   ' 

g      = 

-  Hi 

7 

To  check  how  good  this  is,  we  consider  a  special 


o 


form  for  f .   The  case  f  =  f  je~^^^  ^^  is  studied  in  the 
appendix,  and  this  value  of  V-,  obtained  is  shown  to  be 
a  very  good  bound. 
Our  result  is 

-  VT(f,Vf)  <  -2-  (f,Kf)^/2  , 

V-,    =  T^  V4egk7      . 
Using   equations    (23)^    we   have 


j^  =    yi-ri^  E^  +  T^Hf 


Ey=<H>=    (1-T]^)E   +  T]^(f,Hf)    ,         as    (^,f)=0 

Setting  ■    (f,Hf)    =    (f,Kf)    +   (f,Vf),    we  get 

Ti-^(f,Kf)    -  -^  ^    (f,Kf)^/2    <   Ej.   -    E    +  Ti^E    . 
y2      ^1 
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¥e   set  T]^   =(E   -E)/(E^-E)    to  get 

(f,Kf)    -  -1-  -^    (f,Kf)^/^   <   E, 

y2  ^1  -  -  1 

Solving,    we  get 

(46a)  (f,Kf)^/2   <  a    +    yi^+E^  =  kJ/2 


(46b)  "■ 2_ 


ot 


2y2  V^ 


Kp^  turns  out  to  be  several  times  as  large  as  the 
expectation  value  of  the  kinetic  energy  for  the  polarized 
wave  function. 


6.   A  Variational  Principle  for  Quantities  Other  Than 
the  Energy. 

A.    General  Formulation. 

The  following  variational  principle  has  been  set  up 
to  see  how  good  the  polarized  wave  function  is  for  various 
quantities  of  physical  interest.   Suppose  we  want  to 

study  the  physical  quantity  Z.   Z  may  be  the  mean  square 

ikr 
radius,  the  electric  form  factor  e    ,  the  Coulomb  energy 

e^/r,  the  matrix  element  for   the  beta  decay  of  the  triton, 

the  tensor  force,  etc.   We  write  AZ,  the  difference  between 

the  correct  expectation  value  of  Z  and  the  expectation 

value  calculated  with  the  approximate  wave  function,  i.e., 

the  polarized  wave  function,  in  terms  of  a  Green's  function, 
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A  variational  principle  for  this  Green's  function  Is 
then  derived,  and  is  used  to  estimate  AZ .   ¥e  say  a 
quantity  is  derived  from  a  variational  principle  if  it 
can  be  expressed  in  a  form  involving  an  approximation  to 
the  quantity  of  interest,  which  is  stationary  with  respect 
to  small  changes  in  the  approximate  quantity.   That  is  to 
say,  if  we  make  a  small  change  in  the  parameters  of  the 
approximate  quantity  by  an  amount  e,    the  difference 

between  the  correct  value  and  the  variational  expression 

2 
is  of  order  e  .   This  means  that  the  change  does  not 

depend  upon  the  sign  of  e,    and  so  is  stationary. 

An  example  of  a  variational  principle  is  given  by 
the  Rayleigh-Rltz  variational  principle  for  the  energy. 
A  wave  function  which  minimizes  the  energy  is  a  good 
approximation  to  the  actual  wave  function,  and  <H>  is 
stationary  with  respect  to  sm.all  changes  in  this  wave 
function.   When  we  say  that  the  energy  is  derived  from 
a  variational  principle,  we  mean  that  we  choose  the 
energy  as  that  energy  which  is  derived  from  a  wave 
function  which  makes  <■  H  )>  stationary .   The  Rayleigh-Rltz 
principle  is  also  a  minimal  principle,  in  that  it  tells 
us  how  to  get  the  best  value  of  the  energy.   A  variational 
principle  does  not  have  this  property  in  general .   It 
merely  has  the  statlonarlty  property. 

The  property  of  statlonarlty  is  best  illustrated 
by  the  vertex  of  a  parabola.   As  we  m.ove  slightly  away 
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from  the  vertex,  the  change  In  the  height  of  the  parabola 
Is  of  second  order  In  smallness.   Furthermore,  the  height 
always  changes  in  the  same  direction,  whether  we  move  to 
the  right  or  left  of  the  vertex. 

A  variational  expression  for  a  quantity  is  a  Taylor 
series  expansion  for  the  difference  between  the  correct 
value  and  approximate  value,  truncated  after  the  second 
term.   In  this  case,  the  difference  between  the  function 
and  the  variational  expression  does  not  depend  upon  the 
sign  of  the  expansion  parameter,  and  so  is  stationary. 
It  is  not  necessary  to  have  some  type  of  integral,  as 
one  does  in  the  calculus  of  variations,  in  order  to 
speak  of  a  variational  principle  --  the  only  important 
concept  is  the  concept  of  stationarity .   This  type  of 
variational  pi^inciple  can  be  used  as  an  Iterative  means 
of  calculating  the  function.   We  guess  at  a  zeroth 
approximation,  put  it  in  the  variational  expression, 
evaluate  a  better  approximate  value  of  the  function, 
and  Iterate.   This  is  the  Newton-Raphson  method  for 
solving  an  algebraic  equation.  ' 

In  deriving  a  variational  principle  for  the  wave 
function,  we  note  that  if  we  write  the  true  wave 
function  ip   in  terms  of  the  approximate  wave  function  ^   as 

^  =  (j)  +  A  (  .  .  .  G ) 

and  go  to  first  order  in  A,  we  have  a  variational  principle 
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for  the  difference  between  the  true  wave  function  and  the 
approximate  wave  function.   If  we  define  ^   In  this  means, 
any  expression  Involving  (j)  would  be  variational.   The 
term  Involving  A  Involves  the  Green's  function  G,  which 
can  be  expressed  In  terms  of  a  variational  principle,  as 
shown  below.   Substituting  this  Green's  function  in  AZ, 
we  have  a  variational  principle  for  AZ .   We  then  proceed 
in  a  way  analogously  with  the  Rayleigh-Rltz  variational 
principle.   We  guess  at  some  forms  for  the  unknown 
functions  involved  in  AZ,  and  vary  parameters  until  a 
stationary  point  is  found .   We  then  have  AZ  expressed 
in  terms  of  a  variational  principle.   If  the  resulting 
value  for  AZ  is  small  compared  with  Z  ,  we  can  then 
say  that  our  wave  function  is  quite  probably  a  good 
wave  f\j.nction  for  the  quantity  Z. 

The  variational  principle  is  derived  as  follows. 
Let  ^   be  the  normalized  wave  function  found  to  give  the 
best  value  of  the  energy,  i.e.,  the  polarized  wave 
function  of  Section  p.      Let  this  energy  be  denoted  by  K, . 
We  again  define  expectation  values  with  respect  to  (}), 
that  is, 

<Z>=  ((^,Z<i))/(^,(l)) 

for  any  operator  Z.   Using  this  notation,  E^^  is  defined 
as  E  =<^H'^.   Let  E^   be  a  hamiltonian  of  which  (j)  is  an 
eigenfunctlon; 


(^7)  Hq<1)  =  E^.c}) 
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^  and  E  are  defined  above  in  equations  (l8) .   Using  the 
three  body  correction  potential  of  Section  4,  W,  we  have 

(48a)  (Hq-W+F)^  =  Ey^  , 

(48b)  F  =  E^-E 

We  now  write  f   in  terms  of  (j)  in  first  order  perturbation 
theory.   Since  a  variational  principle  is  to  be  constructed, 
there  Is  no  need  for  further  terms. 

(49)  f  =   ^   +  p,  -^tT   [-W+P-<-W+P>](l)  . 

In  equation  (49),  the  term   <( -W+P_^  Is  subtracted  as  we 
are  Interested  only  in  the  principal  part  of  the  Green's 
function.   In  view  of  the  fact  that  <^W)>  =  0,  this  simplifies 
to 

(50)  ^  =  <t'  -  -w^-  "^'^    ' 

Since  we  neglect  second  order  effects,  we  have  (?//,!//)  =  ((t),(j))=l 
The  quantity  t]  of  the  previous  section  is  a  second  order 
quantity  and  is  neglected. 

We  proceed  to  find  a  variational  principle  for  the 
Green's  function  l/(Ey-HQ)  as  follows.   Let  A  =  E^-Hq . 
Consider  an  operator  B  %  A 


^  =  -TT-^ TT- 

B   +(A-B   ) 


B"-^    B"^         B  -" 
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or 

(51)  G  %  2B  -  BAB  . 

Since  this  is  correct  to  second  order,  it  is  a  variational 
principle o   This  Green's  function  has  the  further  property 
that  it  is  negative  definite,  if  we  assume  that  E^^  is  the 
lowest  energy.   This  makes  it  possible  to  view  equation 
(51)  the  following  way.   Let  A  =  -(Kj-H^)  be  positive 
definite.   We  have 

(A"-^-B)A(A"-^-B)  >  0  . 

We  solve  this  as  follows: 

(A"-'--B)A(A"-'--B)  =  (A"-^-B)(l-AB) 

=  A"-*"  -  B  -  B  +  BAB  . 


We  get 


-  G  =  A  -^  >  2B  -  BAB 


The  quantity  2B-BAB  has  a  fixed  upper  bound.   Consequently, 
small  changes  in  the  quantity  (B-A~~)  will  produce  changes 
in  2B-BAB  independent  of  the  sign  of  the  error  of  (B-A~  ) . 
This  means  the  expression  is  stationary  about  the  correct 
result,  and  so  is  a  variational  principle. 

We  define  the  projection  operator  P^  as  the  projection 
on  the  trial  state 

(52)  Pq  =  \^y<}\  . 

Note  that  this   projection  operator  is  not  the  same  as  the 
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projection  operator  of  equation  (20)  above.   We  likewise 
define  Q^  =  1-Pq,  and  express  W  as  In  equation  (10): 

(55)  -  W  =  PqHQq  +  QqHPq  . 

Using  this  variational  expression  for  the  Green's  function, 
we  get  for  a  variational  principle  for  AZ 

(5^)         AZ  %  -  2  Re  (zQq[2B-B(Eq-Hq)B]QqW/>  . 

It  is  not  necessary  to  know  the  entire  matrix  B  -- 
Instead,  the  function  |B^^  and  {b^^   are  sufficient,  where 
they  are  defined  as  follows.   We  first  note  that 

QqWcJ)  =  Wcf) 

from  equation  (55) ^  noting  that  Qq  =  Qq,  as  Qq  Is  a 
projection  operator.   We  also  note  that 

QqZI^)    =   Z|(t))  -  \^){^\Z\<\>) 
or 

QqZ|(|))  =  (Z  -  <Z>)  Icf))  . 

We  notice  that  the  only  vectors  occuring  in 
equation  (54)  are 

(55a)  IB-^)  =  B(Z  -<Z>)|([)) 

(55b)  Ib^)  =  -  BW|(1))  . 

If  we  multiply  equation  (55a)  by  -((t)jw  we  get 

t56a)  -  ((t)|w|B3_)  =  (B2l(Z  -<Z>)|^)  . 
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Equations  (55)  and  (56a)  define  the  vectors  |b.).   We  can 
Impose  additional  conditions  on  |b.)  as  follows.   We  note 
that  the  Green's  function  occurs  multiplied  by  the  projec- 
tion operator  Qq  to  the  complement  of  the  space  of  |(f)). 
Since  B  is  to  approximate  G,  we  require 

(56b)  ((t),B^)  =  0  ,         i  =  1,2 

Substituting  equations  (55)  into  equation  (5^)  we  get 

AZ  %  2  Re[-(B^,¥(j))  +  (cf),  (Z-(Z^)B2) 
-  (B-^,(Ey-HQ)B2)]  , 

subject  to  the  restrictions  (56).   Now 

-  (B-^M)    =    (B^AqH^) 
=    (B-L,H(t)) 

by  virtue  of  equation  (56b).   Likewise,  (B, ,WBp)  =  0, 
and  so 

(B-l,HqB2)  =  (B^,HB2)  . 

Let 

(57)  %^  =  H  -  Ey  . 
We  get  for  AZ: 

(58)  AZ  ^   2   Re[(B^,H(j))  +  iB^,Z^)    -    {B^,>Ib^)]     . 

We  notice  that  the  B.  occur  linearly  except  for  the 

t 
combination  B-,Bo-   We  therefore  transform  B.  to  B. 

as  follows 
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(59a)      b|  =  |(B-L  +  Bg)  ,  Bg  =  |(B-j_  -  Bg) 


or 


(59b)      B-^  =  BJ  +  Bg  ,   ^2  =  ^1  ~  ^2  '   %^2  =  B^  "  Bg  • 


We  get  for  AZ: 

(60)  AZ  ^  2[(B^,^BJ)  +  (B-[,H(1))  +  (B^^Zcf))] 

-  2[(B2,^B2)  -  (Bg.Hcf))  +  (B2,Z(t))] 
subject  to 
(6la)      FQ  =  (B^,H(t))  -(B^,Z(i))  +  (B2,H(t))  +  (Bg.ZcJ))  =  0 

(6lb)  {b[_A)    =   0    ,  1  =   1,2. 

We  note  that  (B.,aB^)  >  0,  as  E  is  the  lowest 
energy  found  by  the  minimum  principle  for  the  energy. 
This  is  used  to  help  us  find  a  stationary  point,  as  follows 
We  Introduce  a  Legendre  parameter  A,  and  look  for  the 
stationary  point  of  AZ  +  AFQ,  with  A  chosen  so  that  PQ  =0 
at  the  stationary  point.   Set  A,Z  =  AZ  +  APQ.   A,Z  is  the 

sum  of  two  quantities,  which  will  be  called  RMS^,  where 

t 

the  subscript  i  refers  to  the  subscript  of  B. .   Since 

I  RMS.  I  is  bounded  from  below,  a  stationary  point  for 
A  Z  will  be  a  saddle  point  --  a  minimum  for  RMS^  and  a 
maximum  for  RMSp-  See  below  for  details  on  finding  the 
stationary  point. 
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The  negative  deflnitiveness  of  the  Green's  function  plays 
a  significant  role  in  helping  us  reach  a  stationary  point 

B.   Specialization  to  Error  in  the  Mean  Square  Radius. 

2 
The  special  case  Z  =  r  will  now  be  studied.   An 

I        2 
ansatz  is  chosen  for  the  B.,  and  Ar   is  varied  subject 

to  the  restraints  until  a  stationary  point  is  found. 

I  2 

Although  the  choice  for  B.  is  arbiti-ary,  Ar  will  satisfy 

a  variational  principle,  as  first  order  variations  in 

the  parameters  of  B.  will  cause  second  order  variations 

m  Ar  . 

In  choosing  an  ansatz  for  the  B.,  we  notice  that 

they  approximately  satisfy  the  equation  (Erj-Hp^)B.  =  0 

at  large  distances  from  the  origin,  when  ^   is  close  to 

zero.   B,  approximately  satisfies  this  equation  very 

close  to  the  origin,  and  Bp  in  the  region  of  the  mean 

square  radius,  where  W  is  close  to  zero.   Consequently, 

the  ansatz  for  B.  will  be  a  function  like  (|),  with  a 

different  set  of  parameters: 


2         2f,,  '     3   ji- 


^  (x2+4a,,y2) 
(62a)      b!  =  b.  'V~   c  ..  e    ^^ 

+   symmetric  permutations  of  v-.  >_rp,r^; 

(62b)  ^2i  =  ^  • 

The  parameters  c^  .  are  determined  by  the  orthogonality 
relations  (6lb).   The  b,  are  determined  as  follows. 
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RMS.  can  be  expressed  in  the  form 

(63)  RMS^  =  b^  X  FN^  +  b^PM.  , 

where  FN.  and  PM.  are  functions  Independent  of  b..  We  require 

Srms. 

1 
We  obtain 

FN. 

(64)  b. 


1 


i      2X-PM .  ' 

The  procedure  to  find  the  non-linear  parameters  Is  as 
follows.   We  first  determine  the  linear  parameters.   The 
quantity  RMS^-RMSp  Is  then  minimized.   This  will  bring 
us  near  the  stationary  point.    We  then 
fix   .  A  at  some  value,  and  look  for  the  minima  of 
|RMS.  I,  and  then  try  different  values  of  A  until  we  get 
a  result  converging  on  PQ  =  0 . 

In  order  to  do  the  minimization  on  the  computing 
machine.  It  Is  necessary  to  have  some  Idea  of  the  size 

of  the  parameters.   The  limits  on  the  sizes  of  the 

t 
parameters  contained  In  B.  are  guessed  at  by  physical 

reasoning  --  they  probably  will  not  be  very  much  larger 

than  the  parameters  of  ^.      There  remains  the  problem  of 

finding  the  size  of  A.   There  are  two  ways  of  doing  this. 

One  way  is  to  notice  that  RMS.  can  be  written  as,  using 

equation  (64) 

2 

^  PNT 

(65)  RMS^^  =  -  ^  RM^  . 
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and 

This  form  suggests  that  |a|  <_   2.   This  restriction  on  A 
has  helped  considerably  In  finding  a  stationary  point. 

Another  way  to  determine  A  Is  by  doing  so  directly, 
We  write 


{66) 


PQ  =  2  b^  PQ^ 


Using  the  requirement  PQ  =  0,  we  get 


(67) 


bg  =  -  b^PQ^/PQ2 


We  write 

(68) 

Equation  (64)  becomes 

(69) 


PN.  =  FY^  +APQ^ 


bg^  - 


2  FM^ 


Equations  (6?)  and  (69)  are  two  Independent  equations  which 


can  be  solved  for  7\  in  terms  of  b- 


We  use 


(70) 

and get  for  A 
A  =  - 


FY^-AFQ^ 


PMg  PY^       FYg 
2  ^^1  -  FQT 


FM 


1  FQ2    -^    ^"=-2 


1  + 


FMg/  PQ;l  \' 


This  serves  as  a  guide,  telling  us  how  far  we  are  from 
the  stationary  point. 


50   - 


2    2 
The  result  is  Ar  /(v  )>  =   +.OI57,  showing  that  to 

a  very  good  probability,  the  polarized  gausslan  wave 

function  Is  a  good  wave  function  for  the  mean  square 

radius  as  well  as  the  energy. 

7.   Tables  of  Results  . 

A .    Triton 

We  have  the  following  results  for  the  minimum  in 
the  energy,  calculated  with  a  polarized  gausslan  wave 
function.   The  potential  is 


v  =  ^0^ 


2 

-Ltr.  . 

e    "^   mev. 


1<J 


2  2 

The  value  of  the  constant  -fi  /m  =   41.5897  mev  f  except 

where  mentioned  otherwise.   We  define  the  following 

symbols 

E  =  <^H  ^  in  mev. 

R  =  <^r'^y.-^'^'^   in  fermi 

Let  Y  be  the  Coulomb  energy  (for  He^) .   In  studying 
the  Coulomb  energy,  the  value  of  m  was  not  changed,  so 
that  the  effects  of  the  Coulomb  force  are  more  readily 
seen.   The  param.eters  of  V  are  Identified  as  follows. 

Laskar^'   takes  V^  =  -37-44  mev,  [i  =    .2669  f~  • 

n  2 

Tang  et  al .  take  Vq  =  -51-5,  1^  =  1/(1. 6f)  . 

Baker  et  al.  have  taken  this  same  value  of  [x,    and 
Vq  =  -48.0,  -51-5,  -55-5  mev. 
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For  Laskar's  potential,  we  have  the  following  results. 
Using  the  best  polarized  wave  function,  we  get 

E  =  -8.496  R  =  1.65 

Introducing  a  Coulomb  force,  we  get 

E  =  -7.806       R  =  1.68       Y  =  .684 
For  the  unpolarized  gaussian,  we  get,  with  Y  =  0 

E  =  -6.675  R  =  1-^8 

With  one  polarized  wave  function,  we  get,  where  the 
variational  parameters  are  f  and  b  =  4a/3: 

f  =  .618        b  =  .430       E  =  -7.57 
f  =  1.69        b  =  3-^3      E  =  -7-51 

2  2 

Using  Tang's  value  of  -fi  /m  -  41.47  mev  f  ,  and  the 

best  polarized  wave  function,  we  get 

E  =  -8.571  R  =  1.65  . 

The  results  for  Tang's  potential,  we  have: 
With  the  best  polarized  wave  function 

E  =  -9.52  R  =  1.46 

For  the  unpolarized  gaussian 

E  =  -6.80  R  =  1.285 

E  =   0.84  R  =  4.66  . 

For  a  gaussian  wave  function,  there  are  two  values  of  the 
variational  parameter  which  yield  minima  for  <H_^.   The 
second  minimum,  however,  is  not  a  bound  state. 
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For  a  sum  of  two  unpolarlzed  gausslans 

E  =  -8.98  R  =  1.41 

For  a   single  polarized  wave   function: 

f  =    .618  b   =    .430  E  =   -7.40 

f  =   1.69  b  =   3.^3  E  =   -6.96 

We  see  that  most  of  the  improvement  is  caused  by  the 
resonance  effect,  that  is,  by  taking  a  sum  of  functions. 
However,  the  polarization  effect  is  quite  large.   The 
effects  of  resonance  and  polarization  are  about  of 
equal  importance  in  evaluating  the  energy  or  the  mean 
square  radius.   For  Tang's  value  of -fi  /m  and  for 
the  best  polarized  wavefunction,  we  get 

E  =  -9.62  .  R  =  1.45  . 

For  Baker's  potentials  we  get 

V  =  -48.0   E  =  -6.71   R  =  1-58  Eg  =  -6.78  Eq  =-4. 06 

V  =  -51.5   E  =  -9.52   R  =  1.46  Eg  =  -9.42  E^  =-6.9^ 

V  =  -55.5   E  =  -13.09  R  =  1.35  Eg  =  -12.15  E^  =-9.64 

Et-,  is  the  energy  which  Baker  et  al .  obtained. 

D 

Eg  is  the  best  energy  obtained  with  a  gaussian  wave  function. 

B.  jpper  Bound  to  the  Mean  Square  Radius. 

We  work  with  the  potential  of  Tang  et  al .   Define 
the  dispersion  in  the  energy 

HH2  =  [<h2>  -  <H>^]  /  (E-^  -  <H»  -  <H>   . 

For  the  best  polarized  wavefunction,  setting  E  =  -lO.mev, 
we  get 
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<H^>  =   107-59  xnev^       t\^  =  .05  V^  =    .182 

Kq  =   13^.3  mev  HH2  =  11. 58  mev        Gq  =    .11  mev"-"- 

<'Wr^V{>=   57.21    (mev  f)^  <(  w^J>  =   17-17 
Ar^  =   10.28   f^ 


2 

mev 


For  the  best  unpolarlzed  wave  function,  we  get,  with 
E  =  -10 .  mev: 

<H^>  =  260.76  mev^    T\^  =    .333    HH2  =  40.^0  mev. 
Gq  =  .16  mev" 

<Wr^W>  =   433.01  mev^f^        <¥^>  =   224.70  mev^ 
Ar^  =   114.42   f^ 

For   the  best  polarized  wavefunction,    but  E  =   -9-521,    i.e., 

Ey   -   E  =   10"^    ,      we  get     Ar^  =   3-184 

I   2 1 
Fixing  E  at  this  value  and  minimizing  |Ar  |,  we  get 

Ar^  =  3.16. 

The  upper  bound  for  the  error  in  the  mean  square 

radius  is  of  order  several  times  the  mean  square  radius. 

Hence  this  upper  bound  is  of  no  help  here.   It  would  be 

useful  in  situations  where  the  dispersion  is  very  small. 


C .  Variational  Principle  for  the  Error  in  the  Mean 
Square  Radius. 

The  result  for  Laskar's  potential  by  using  the 

variational  principle  is   p 

-^^  =  +  .01573. 
<r  > 

B  represents  an  extremely  compressed  state,  as  the 

quantity  a.  .  is  around  8. 

The  program  was  run  in  double  precision,  as  some 
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quantities  were  very  large.   For  example,  FM  =  -2  x  10  . 
The  results  are  (refer  to  equations  (62)  for  the 

definitions  of  the  parameters): 
^ll=~^-°^5'  ^1=  -^-96^  X  10"^,  c^2=-15.0678,  b2=-. 2936x10"^ 
f^^  =  0.8437,  a^-j_  =  7.245,  f^^  =  1.007,  ag^^  =  5-086 
f^2  =  0.6022,  a^2  =  9-087,  f^^   =  0.697,  a^^  =   2.255 
A  =  -.0310,   FQ  =  -.45  X  10"5. 

2 

A  slight  change  in  A  alters  the  parameters,  but  not  Ar  , 

as  this  is  a  stationary  point. 


D.    Results  for  Alpha  Particle. 

With  Laskar's  potential,  we  have  the  following  results 
For  a  non-polarized  gausslan,  we  have 

R  =  1.53     E  =  -29.36     7  =  .3265  f~^ 

For  a  single  polarized  gausslan,  we  have 

f  =  1.94,  -^  a  =  2.27,  ^  b  =  2.96  (compressed  polarization) 

R  =  1.33  E  =  -30.16 

f  =  0.496,-^  a  =  0.338,^  b  =  0.329  (extended  polarization) 

R  =  1.3^5  E  =  -30.49 

f  =  1.25,  -^  a  =  2.49,  -I  h  =  0.711  (mixed  polarization) 

R  =  1.3^5  E  =  -30.45 

For  a  linear  combination  of  the  compressed  and  extended 
polarizations,  we  have 
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3 


c^  =1  f^  =   0.412       ^  a^   =   0.272         "I  b^   =  0.567 

Cg  =.106  fg  =   2.955       1^2"   2.266         I  1^2  "   2.117 

R  =   1.580        E  =   -51.55       Y  =    .885 


For  a  linear  combination  of  all  three  types  of 
polarization,  we  have  E  =  -51.40. 

Kalos   has  used  the  following  parameters  for  the 
potential:   Vq  =  -59-6,  [x  =  .512.   He  has  set  the  Coulomb 
potential  equal  to  zero.   He  obtained  an  energy  of  E  =  -28.5 
mev.   With  a  polarized  wave  function  of  7  parameters,  the 
following  is  obtained: 

R  =  1.56       E  =:  -27-84  Y  =  0 

With  an  11  parameter  function  : 

R  =  1.57       E  =  -27.95 
The  improvement  in  the  binding  energy  obtained  by  using 
a  wavefunction  which  is  a  linear  combination  of  the  three 
types  of  polarizations  is  only  .09  mev  better  than  the 
wave  function  which  is  a  sum  of  the  compressed  and 
extended  polarizations  only. 
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APPENDIX 

1 .   Physical  Constants. 

—27 
■fi  =   1.05450  X  10      '    erg   sec.      See  Reference   I5. 

1   amu     =   1.65979  X  10"^     gm 

1  mev  =  1.60210  X  10"°  erg 

14   ^ 
To  get  m,  we  use  the  tables  in  Lelghton.    H-^  has  a  mass 

defect  of  [I5.835  +  .005]/931 .l4l  mev.   The  electron  mass 

is  .510976/931.141  mev.   He^  has  a  mass  defect  of  I5.817 

amu  ,  and  He   has  3 -607.   We  subtract  the  electronic 

mass.   For  H-^  we  have  for  m 

m   =  [1  +^  (15.835-.511)/931-I4l]  amu 

m   =  1  + -^  (15.817-  1.022)/931.l4l 
He^      ^ 

m  4=  1  + -F  (3.607  -  1.022)/931-1^1  . 
He         - 

13 

The  charge  of  the  electron  is  -^ 

e   =  4.80298  X  10"-^°  statcoul".  . 


The  Triton  Problem. 


We  define  our  coordinates  x,  y,  R  as 


(Al)      2S  =  Zi  "  £2 

(A2)     y  =  n^  -  |(i:i  +  L2^ 

(A3)      R  =  3  2  r^ 

The  following  relations    hold   between  r^,  x  and  y: 
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(a4)  r^   -  1:2  "  —  ~  -^ 

(A5)  Zi   -  1:3  =  ^x-Z  =  2S2 

(a6)  £2  "  -3  "  ~\^'  Z  =  2Si 

(AT)  £3  -  Idii  +£2^   =  Z  =  I3 

(A8)  £2   -   ""^  g  "^  =  ^  [^2   "  -5  "  -1   +  -2^    "    "  ^  ^'  i^  "  ^2 

^2   +  ^^        ^  1 

(A9)  r^ 2 — ^  =^2S--^Z  =  Zi    • 

¥e  can  solve  for  r.  as  follows. 
Z  -  5R  =  -  I  (Zi  +  r^) 
I  (Z  -  3R)  +  X  =  -  Srg 

(AlO)    i:2"~i-"3^^- 

(All)     r^=|x-3Z+R 

Z3  =  5R  -  Zi  -  1:2 
(A12)     Z 3  =  R  +  I  Z 

Define  £  .  =  r.  -  r . 

7-  r?  .  =  5Z  X? 

IZ  Z- J  =  I  x2  +  2y2  . 

2 
The  square  radius  r   is  defined  as 

(AI3)     r^  =  i^  (r^  -  R)2 
This  can  be  expressed  as 
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r'  =  ^(-ix-iz)2  4<-i^-32''4'|i'' 


or 


(A14)  r^  =   1    (x2  +1  y2) 

We  also   have: 

(A15)  ^Zl  Zij  =   9r^    . 

We  write   our  polarized  wave   function   as 

^^^  =   exp    [-  -^   (xj   +  a^jj)]    . 

If  we  set  f.  =1,  a.  =  4/3,  7  =  -^  A,  we  get  Laskar's 
trial  wave  function,  with  A  Laskar's  variational  parameter 

(t)-L^  =  exp  {  -  2^  [(-^  +  y)^  +  a^(|  X  +1  Z)^]}  • 
2  i 

The  spatial  part  of  the  wavefunctlon  of  the  trlton 

can  be  written  as  cj)  times  a  function  of  R.   We  set  this 

function  of  R  equal  to  unity.   Consequently,  we  write 

for  the  kinetic  energy  operator  T,  where  k  =  -ti  /2m: 

^   2 

T  =  -  k  2^  Vf  ] 

1=1  -  ^  R=0 
2 


V  =  Vq  ^Z  e 


-|j,r 


ij 


Since   (]),    T,    and  V  are   syTTjnetric    in  r^,    we  relabel 
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the  indices  in  the  integrations  of  {^,Q^),   where  9  is 
any  operator,  so  that  only  cj)   appears  on  the  right. 
Furthermore,  since  (j)  .  is  invariant  under  the  transforma- 
tion rg  -^  Li'    we  also  have  (^^.  ,Q^      )   =    {<\>^.  ,Q^      )  . 
Consequently,  the  evaluation  of  {^,Q^)    reduces  to  the 
evaluation  of 

An  expression  for  Tc()   will  now  be  derived. 

where  r  ..  is  the  i-th  component  of  the  position  vector 
of  the  j-th  particle,  r.  .  is  not  to  be  confused  with  r.  .. 
Let  X.  be  the  i-th  component  of  x.   This  is  not  to  be 
confused  with  the  vector  x. .   We  have  the  following 
relations,  where  S.  .  =  S/Sr.  .: 


^li  ^±  =   ^    '  ^21^1  "  "^  '  ^^1^   i  =  0  , 

^li^i   =  -  i  '     ^21^1   =  -  i  '      ^31^1   =  ^  • 


ky 


Tcl)^^  =   #:2Z  [  S^j.(2x-a,y)  .  +  S^  .  (-2x-a.y)  . 


+  ajj  2a^yj]  e 


5k  (-'+-1^") 
e 


(A16)   T*,,  =  ^[^,-1^-2  ^/_|  Z^  a^y^]*5. 
The  integrals  involved  are  evaluated  below. 
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5 .   Minima  in  ^H/»  for  a  Gaussian  Wave  Function . 

The  expectation  value  of  <H)>  for  a  gausslan  wave 
function  Is  a  simple  expression,  making  It  easy  to  study 
values  of  the  parameter  of  the  wave  function  which 
minimizes  (U).      We  work  with  the  wave  function 

T^  =  I2k.j{l-yv^)(\>    . 

To  do  an  Integral  like  (([>,({)),  replace  the  volume  element 
dx  dy  by  dr.   We  get 

3 


i^A)    =Je-^y^      6r=    (^ 


r^  =  TrV^^^^  6r/{^A) 


1  A  (^,^)/i^,^) 


(A17)  <^!>  =  ^ 

<T>=  12k7(l-(7r^>) 

(Al8)         <(t>-  6k7  . 

To  get  <(V/»,  we  note  that  <(v)>   can  be  written,  due 
to  symmetry,  as 

V   =  5Voy'e-^^^^-^^^W(<l>,ct)) 

=  3Vq  /  e-^^^  -^^^^^y  ^dx  dy/Je      ^  ^^^^ 


-   61 


Now  ,  P  ii  p 


^   dx  dj  =  (5  /|)^ 


7 

And  so 


e 


^   dx  dy  =  Tr^(^)^/^  (7-Hi.)  ^^   7  ^^ 


We  therefore  get: 


or 


Let 


We  get 


(A19)  <V>=   3Vq(1  +^)"^/^ 

(A20)  <H>=    6k7    +   3Vq(1    +^)"^/^ 


dH  .   6k  +  3V^|(1   +it)-5/2ii_ 
d7  ^0    2^  7  2 


fi=  0    If 

d7 


61c/  =-  Vol     (l+tt)-5/2 


,V,2 
a    =    Itl — ) 

3*^ 


(A21)  (7-h-L)"   =  ay 

Set  f(7)  =  (7-I1-L)-  -  ay.      To  study  the  zeros  of  f(7),  we 
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look  at  the  derivative  of  fiy) 


f'(7)  =  5(7Mi)^   -  a 
Set  f ' (7^)  =0.   We  have 

72  =  (3  ex)     -  M-^   as  a  >  0,  |j,  >  0  , 

and  we  are  Interested  only  In  positive  values  of  7. 
For  the  potential  of  Baker  et  al .  we  have 

M-  =  .3906,     a  =  .5264,     7^  =  .178  . 

¥e  know  one  zero  of  f(7),  which  will  be  called  j^ 

f(7o)  =  0 

7q  =  .301  ,         E  =   -6.803  mev  . 

Since  f(0)  =  ij,-'  >  0  and  f'(7„)  >  0,  there  must  exist 
another  zero  of  f{y^),    such  that 

0   <  -y^   <  y^    . 

Since  f'(7)  >  0  for  7  >  y^,    f{y)    >   0  also. 

We  see  that  there  are  only  two  zeros  of  f (7) .   The 

result  Is  y.    =    .023,  E  =  +  .844  mev. 

This  shows  that  If  there  Is  an  excited  bound  state, 
It  Is  very  poorly  represented  by  a  gausslan  wave  function. 

The  above  zeros  were  found  by  the  minimization  routine 
SC00P . 
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4.  The  Error  In  the  Energy  Due  to  Neglect  of 
Relatlvlstlc  Effects. 

The  order  of  magnitude  of  relatlvlstlc  effects  will 
be  calculated  here. 

The  expectation  value  of  the  kinetic  energy  of  the 

trlton  Is  about  25  mev.   This  means  each  particle  has     ^ 

an  average  energy  of  8.33  mev.   We  calculate  7  =  (1-  -^)   ^^ 

c^ 

Since  the  kinetic  energy  Is  proportional  to  the  mass,  we 
can  assume  that  relatlvlstlc  corrections  to  the  energy 
will  be  about  the  order  (7-I)  times  the  energy. 

vf  _  2  X  8.^3 

2  "      2 
c       mc 

=  .0178 

(7-1)  =  .00886 

(7-I)  times  25  mev  =  .666   mev. 

2 
The  relatlvlstlc  corrections  are  of  order  ^  mev. 

5 

5.  The  Vlrial  Theorem. 

The  vlrial  theorem  is  as  follows .   We  define 


(A22)  rV  =  >   L±    '  1±    ' 

The  vlrial  theorem,  states 

(A25)  <Ct  -  I  r[V,V])>=  0  . 

To  test  this,  we  need  the  quantity  -^  rVV.   We  have 
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2 
V  =  Vq  >   e    ^ .   We  see  that  V  is  the  sum  of  terms  like 

^12^^1'^2'^3^  =  Sl2''^12^"   '^'^'^^  ^^' 

2 

-ur.  . 

We  have  the  following  relations: 


^^12        ^^12  ^^^2  ^2^2  ^^12 

"5rZ~  "   dr-j^2    '  "SF^  "    ~   dr.  g    '  8r^ 


rVf-L2   =    ^^1  •  ^12    -  ^2  •  ^12^g  =  ^12  "  ^1 


2  '^"'^12 

-12  '-12^   "    "^-^^12    ^ 


=    0 


(A24)  -i  rVV  =   -^lVq  ZI  ^i^-   ©        ^"^ 


If  V  includes   a  Coulomb    energy  V",    we   have 


^12  •  ^12^12  =    "    Si2    • 


(A25)  I  rVV"^   =    -  "I  V^ 


The  following  theorem  is  proven- 
Theorem.   If  a  variational  function  permits 
a  scale  transformation,  the  value  of  the  variational 
parameter  which  yields  a  m.inimjam  to  the  energy  also 
satisfies  the  virial  theorem  exactly. 

Proof.   By  hypothesis,  ^'  =  ■'p[l\v-^,'hv^,  .  •  .) ,    is 
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also  a  variational   function.      Define   expectation  values 
with  respect   to   i//',    i.e.,   <9>  =    (^ '  ,0^)/(^ ',?//')  .      We 
want   to   show   that   if  -tt-  <H>  =  0,    then   equation    (A23)    is 
satisfied.      For   convenience,    we   drop   the  prime   from  Tp . 


We   have 


dA       '^ —  ■JAtT        dA        ' 


H  ^  T  -'* 


We   evaluate  -rz-  ^h)> 


dA 


d 


A-^<H>=    [(rV^,H^)    +    (tlJ,ErS/f)]/{i;,f) 


<^ 


{ij,f)'^ 


o   [(rV^,^)    +   (^^rV'^)] 


Now    {ryf,f)    =   -iffS/TTp)  .      We  get 

d<H> 
A  —^ —     =  <[-VrH   +  HrV   +  Vr<H>   -  <H>rV]> 

Vr  =  rV   +  3n,      where   n   is   the  number   of  particles 
We  get 

A  ~<H>  =  <[H,rV]> 
or 

(A26)  ^<H>  =  <[H,    A.]>    . 

Setting  -^  <H_^  =  0,   we  obtain 
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<[H,rV]>  =   0 


Tr  V   =    -    k  >       V7   r . •   V  . 


2  ^  ^,    ,^—      9 


-   k  y^  r     V^  V        -   2k  '2ZL  vf 


or 


Tr  V   =   rTV    +  2T 

<[H,rV]>  =  <    2T   +[v,rV]> 
This    is   equal    to   zero,  proving   the   theorem. 

/-  3      4 

D.   The  Coulomb  Energy  for  He  and  He  . 

The  Coulomb  energy  E  ,  evaluated  with  a  gausslan 


wave  function,  Is: 


~     2 
e^r  e"^^  X  dx 


E   =   ° 


C     CO      2 

r  e-^^   x2  dx 


0 
E   =  2e2  (2)1/2  _ 

C  IT 

For  a  polarized  wave   function,    we   approximate  E^  by 

(A27)  E^=   2e2    (2)1/^  ^non-pol.    ^ 

^  ^  pol. 

where  r      ^   is  the  root  mean  square  radius  for  the 
non-pol . 

unpolarlzed  function,  and  r^^.   that  for  the  polarized 

p  O  J_  * 

wave  function. 

A  value  close  to  the  experimental  result  is  obtained 
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7-   The  Three  Body  Correction  Potential  Written  in 
Configuration  Space. 

The  three  body  correction  potential  W  can  be 
expressed  in  configuration  space  as 

W  .1E-H)i  _ 
9 

The  coordinate  expression  for  W  will  be  evaluated  here. 

I^z|)i.  E  -  V  -^ 

Tcj)^.  has  been  evaluated  in  equation  (Al6)  .   To  get  Tcj) . . , 

2       2 

we  merely  redefine  x  and  y  occurring  in  equation  (Al6) 

2  ^    2    ,-     ^2  ,   ,3    ,  1   ^2 
by  using  equations  (A5)  and  (A8) . 

Xg  +  ay2  =  x^(^  +  ll  "^  +  X'y(-1+  I  a)  +  y^(l+  -^  a) 

We  set  a  =  -^  a^.   Set  ^11  ^j±   =  ^i  • 

^'^i    ,6   ^9!:i^l   _27   r,2^5   2  2. 


+  [x2(l  +1^  -I  a^)  +  x-y(-l+3|  a^) +y2(i+  ^a^)  ]  .t^^^ 


H-  [x2(l  +15  a^)  -  x.y(-l+3^  a2)+y2(i+^  a2)]^^.} 

-^Lr2    2c.T(l). 

(A29)         W  =  E  -  V  ^  e    ^J ^^ — i  . 

KJ  2c^^^ 

W  has  the  form  of  a  sum  of  gausslan  functions  divided  by 

2 
a  siAm  of  gausslans.   For  large  r,   W  goes  to  infinity  like  r 

due  to  the  kinetic  energy  terms. 
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8.   Evaluation  of<H^>. 


t2  .  ^t^\ 


v^  =  Z  V?  . 

—1 


Let  g.  =  2^  ,  k.  =  I  g.  , 

1 

v\,=  O^,  ....  xa^^.  ....  ,e"^"'- 3 '^'', 

where  S.  .  s  ^^ 

ij   dr.  . 


Modifying  equation  (Al6)  we  have: 

2   4.2 
2         "SX  -  3  ky 

U Toz-^j^i.-^  _L  p/^2^2_^  4  i,2_2> 


2 
-gx   - 

e 

4    ,     2 
3ky 

-      --L^l.B-1-lV; 

1"  ^ife 

-?^    -r  ^    K. 

'b^x  =    1 

^2^  =    -1 

ax  =   0 

\ 

s^y  =  -  1 

^2^  =    -  1 

V  =  ^    ( 

/ 

a^x^  =  2 

a^x^  =  2 

:n2   2        „ 
d^x     =   0 

^    for   each 
Cartesian 

:s2   2        1 
\y     =2 

^2   2        1 
^2^      =2 

:^2   2        „ 
a^y     =   2    J 

c 

omponent 

Dropping   the   subscripts   from  g  and   k,    we   have: 
4, 


V 


-r-^  =   [-12(g+k)+8(gV+^  hV)]^  +  8-2-5-2g^ 
^31  -^ 


+  8.  I  •  3-3k^  +  2{ (8.2g^x  -   8-  ^-k^^) 


3 


4  ,  2 


•(-2gx   +'^  kjr)    +    (-8-2g''x   -   8-   yk  j)-(2gx+^  1^) 

+    (8-   |.2k2^).(-2.   I  ky)} 
=  l6[-5(g+k)+2(gV+^hV)]^+l6-6(g^+k2)-8-l6[g\2+^kV] 
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v^<t> 


3-^  =   I5(g^+k2)    +  I8gk  -   g2x2(20g+12k) 


^^ 


-    I6k2y2(|   1,   ^  g)    ^  ,gV    +  f  gVx2y2   ^-  ^   kV 


The  Integrals  are  evaluated  below, 


9.   Evaluation  of  <^Hr^HX 

<Hr^H>=  CTr^T)>  +  2^r^VT)>  +  <r^V^>  ,     as  [r^,V]  =0 

We  want  [r  ,T] 

T  =  -  k  2  V?  .  ,  where  V .  .  s  -^^ — 
ij  ij    dr 


ij 


2   1 


^  =  3^(Zi-R) 


=  1^  fe^  -i^'-i-jl 


[r?.,T]  =  2kr.  .V .  .  +  2V .  .r .  .k 


=  4kr.  .V.  .  +  2k 


[r?,T]  =  4kr.  -V.  +  6k 
—1   ■'     —1  —1 


>1»Z2'T]  =  2k(r^-V2  +  r^'V^),   as  [r^-V  ]  =0,   i  7^  j 


After  we  perform  the  commutations,  we  set  R  =  0.   We  get 
terms  like  Z\'^2   "*"  — 3'— 2'  ^^1^^^  ^^  set  equal  to  -r-^-V^. 
We  get: 


-   70   - 


[r^T]    .^k[4  2r..V..   +6-5   +  2  Zr.  .V.  .] 
=  4  k  2r.  .V.  .    +  4k    . 


We  eet 


[r^,T]    =  4k[^  +  1) 


(A50)     <Hr^H>  =  <r^(T^+2VT+V^)>  -  |  k  <rVT>  -   4k<T>  . 

Note   that    the    first   term  is   not  <r  H  >,      as    [V,T]    7^  0. 
To   evaluate   -  ^  ^  <^rVT>'. 

"""^^^  -       ^fx2  +  a  y2) 


^^^'i'^i         .r    ^2,    2._       2.r,    ^      9.    ^,    27,    2^  3   .2    2 


■k^ 


11  i 


27^   2,    2    ,    5   ^2   2> 
^^(x     +^  a.y    ) 

1 


^31  1 


-jT^^x  +a^y    j  '.X  -t-  T^ 
i 


10.  Various  Calculations  with  a  Gaussian  Wavefunctlon. 

To  check  the  work  and  the  programming,  the  calcula- 
tions were  redone  with  a  gaussian,  i.e.,  a  non-polarized 
wave  function. 
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Let  (j)  =  e"^^-^  .   From  (AI7)  we  have  <:^r  j>  =  I/27 


Likewise,  we  get 


4v    1    d^ 


<r^>  =  4  ^  {^M/{^,^) 
■^       dy 


or 


(A52)  <^>--^ 

Also, 


<;t>=   12k7[l-  <7r^>] 


(A34)  <T>  =   6k7 


7 
(A35)  <Tr^T>  =   (if2a!   . 

¥e   can  also   evaluate  <^Tr  T>  as   follows. 
<Tr^T>  =   <r^T^  +<[T,r^],T> 

(A56) 

<[T,r^],T>  =    -  I   k  <rVT>  -    4k<,T> 
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rVTcf)  =   -   24-3  k7^(l-7r^)r^(t)   -   12   k7^(rVr^)(|) 


as    follows: 


2  2 

rVr     =   2  r      , 


2         1  ? 

rVr     =  4  rV  z(r.-R) 
:>  — 1  — 


=  I  Z   r    (r^-R) -VCr^-R) 


Now  rVR  =  R  =   0.      We   therefore   have 


2        2 
rVr     -  ^     Zr  r .  -Vr  . 
5  -1     -1 

2^2 
=  3Z   r. 

(A57)  rVr^  =    2r^    . 


We   then  get 


or 


rVTcj)   =    -    24   k7^[    3r^   -    37r      +  r^] 


2   2r 


(A58)  rVT(t)  =   -   24   k7  r   [4-37r    ](1) 


P    2 
We   now  want  to   evaluate  <^r  T  )>.      We   have 


.2,  n.^   1  2      ^   „2    ,,    _2 


T   (()  =    -   12   k  7  Z  V^    (l-7r    )(t) 

=    -    12   k^7  Z    {(l-7r^)V^(()   +  2[V^(l-7r^)]V^({) 
+  (t)V^(l-7r^)} 
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2  2 
In  evaluating  V  r  ,  we  must  not  set  R  =  0  until 

after  the  differentiation. 

(A39)  V^.r2=|(r^.-R.) 

(a40)  V^r^  =  4 


V?^  =   -   37(V^r2)(t)   -    (37Vj_r2)^(l) 


(A41)  V^tj)  =    -    127^   +  \2.yv^^, 


as   it   should, 


S    [V^(l-7r2)]V^^  =   -   37  2    [V^(l-7r2)](V^r2)(t) 

2[V^(l-7r2)]Vj_<j)  =   +  372<i)  2  II  r^ 

=   47^r2(t)    . 
We  get 

^=   -12k27[(l-7r2) 127 ( -1+7^2 )+872r2-47]    . 

(A42)  T^tf)  =   12k27^<{5[l6-327r2+1272r^]     . 

Comparing  with  results  for  the  polarized  gaussian,  we  can 

see  that  the  results  are  in  agreement. 
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11.   Some  Integrals  for  the  Triton. 

2  2 

Let  p    -d-,x  -dpX°y-d^y 

{a43)  1  =       e     ^         "^  ^         dx  dy  . 


To  Integrate  I,  we  complete  the  square  of  the  exponent., 
let 

z  =  V^  X  + 


^^ 


and 


2t^ 


We  have 


d     ^' 


2  2    2     2 

d^x  +  dgX-y  +  d^y  =  z  +  ay 


We  therefore  get  for  I: 

■  =  (4)'" 

(A44)  I  =  Tr58(d^d^-d^)"^/2  _ 

We  want  to  do  the  following  integral 

2  2 

P  -d-,x  -dpX-y-d^y     •  p      o 

J=/e^    ^--3   (  A3_x^  +  Agy^ 
+  E^x  +  EgX^y^  +  E^y  +  E|^(x-y) 

+  F^x^  +  F2X^y^  +  F^xV  +  F4y^  +  P3x'^(x-i)^+Fgy^(x.y)2 

,9      (^p      h  h  26     8 

+  G^x"  +  GgX  y^  +  G^xV  +  G4X  y  +  G^y 

+  Ggx^(x-^)^  +  G^xV(x-Z)^  +  ^8^^^-'^^^ 

+  G^(x-y)^}  ^2^  ^y   • 
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Let  B  s  4d-,d-,  -  do 
J-  3    2 


H  =  B-5/2 


A  -   ^ 


c  =  B-1/2, 


To  find  a  term  of  J,    we  differentiate  the  integral  I. 
For  example,  to  get  the  term  F^-,  we  differentiate  I  as 

\hl   I  • 

We  therefore  need  the  following  derivatives: 

a.H  =  6d^c5 
1      J 

a?H  =  60d^c^ 
1       5 

sJh  =  60-l4  dj  c^ 
1  P 

h^E   =  60.14-18  d^  c^-^ 
1  5 

S^S.H  =  60  d,d^  c^  -  6  c^ 
pi        1  P 

S_S?H  =  60 -14  d. d^  c^  -  120  d^c^ 
5  1  13  5 

S^a?H  =  60-14-18-  d,  dj  c-'--'-  -  60-14-3  dj  c^ 
3  1  13  3 

afs?H  =  -  60'l4-3  d,  d^  c^  +  120  c^  +  60  •14-18 •dFdic-'--'- 
3  1  -^3  13 

SgH  =  -  3  dg  c^ 

^2^  =  3-5  dg  c"^  +  3  c5 
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^2^  =  5'5'7-9  dgc-^^  +  5-5-7-3-2  d^c^  +  3-5-3  c"^ 


Now   Bn  c     =    2nd-,c         .      ¥e  use   this   to   evaluate   St^o 
13  1   ci 

a.S^H  -  3-5-14  d^d-,c^  +  3-10  d-,c^ 
Id  d  3  3 

a^SgH  =  3-5-14.18  d^d^c^^  +  3.10-14  £c^ 
a^S^S^H  =   3-5. 14-18  dgd^d   c-'"-'"   -   3-14-10   d^d^c^ 
+  15-12   c"^ 


We  define  the  following  quantities: 


D^  =  A^d^  +  Agd^  -  Eg  +  I  E^ 


D2  =  E^d^  +  E^d^  +  Egd^d^  +^   d^E^  -  2(P2d^+F^d-^) 


2  .  „  .2  ,  ^  _,  ^   ,1  ^2^^ 

+  |(Pgd^+F^d^)  +  2G^  +  3G^  +  I  G^ 

D^  =  F-^d^  +  F^^l   +   d^d^fPgd^+F^d^)  +^  d^(F^d^+Fgd^) 

-  ^CGgd^  +  G^d^  +  d^d^G^) 

+  ^(Ggd^  -  G^d^d^  +  Ggd^  +  I  d^G^) 

o 
D^   =   G^d^  +  G^dJ  +  -^    (Ggd^  +  G^d^d^  +  Ggd^  +  ^  G^d^) 

+  d^d^(G^d^d^  +  Ggd^  +  G^d-^^)  . 

The  result  for  J  is 

(A45)      J  =  17^  8-6   c5[D^+10c%+10-l4c^D^+10-l4-l8c\] 
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12.   Evaluation  of  <^WGr^GW>. 

Equations  (22)  state  that  G  =  QG.Q.  To  evaluate 

the  above  expression.  It  Is  necessary  to  commute  Q  with 

G-j^.   Define  g'  =  (e'-H)""^,  where  e'  ^  E.  We  know  that 
[Q,G-,]  =  0,  as  •!//  Is  an  elgenfunctlon  of  H.   Since  G^ 

occurs  multiplied  by  Q,  we  can  go  to  the  limit  e'  ->  E, 
and  get  the  result 


(A46)  [Q,G^]  =  0  . 

We  use  the  following  commutation  relation: 
(A47)  [A,B"-^]  =  b"-^[B,A]E"-^  . 

r   2    T 

We  need  to  know  the  commutator  Lp  ,r\     .      Now 

2    1  2 

r  =  ^  2  (r--R)  •   This  can  be  expressed  as 

(A48)  ^2  ^  2  ^  ^  ^2  _  1  ^,  ^^.^^_^  _ 

2 
For  brevity  In  notation,  we  define  r  to  satisfy  r  =  rr. 

We  likewise  define  the  momentum. 
(A49a)  P  -  3  2  p^ 

(A49b)  p^  =  2  (p^-P)^ 

We  get 

(A50)  P^  =  2  P^  -  5P^ 

tPln'^k^^  =  -  2i^Pli^lk^-^ 
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[P,r]r  =   -    lt5  |-  3   [    2  r^   -  2  r^] 


or 


(A51)  [P,r]r  =   0    . 

We   therefore   have 

(A52)  [p^,r]r  =  Z    [p^,r]r    . 

Define  r .  .   as 

with  ijk,i  cyclic.   Define  p,  .  likewise.   We  then  have 
(A54)  [p?.,r]r  =  -  |  •2ifTp.  .r.'  .  . 

We  use  the  following  commutation  relations: 
(A55)  [AB,C]  =  A[B,C]  +  [A,C]B 


(A56)       [ABC,D]  =  AB[C,D]  +  A[B,D]C  +  [A,D]BC  . 

We  have  by  using  equation  (A47): 

[G^r]  =  G^[r,E-H]G^  , 


G-,    p 


We  get 


(A57)         [G;L,r]r  =    -  ^  ^  G^  Z   p^j  G^  r^j 

n 

(A58)  [r'   G3_]  =^21lfp:.  G,  . 
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Now  [QG-^Q,r]  =  QG^[Q,r]  +  Q[G,r]Q  +  [Q,r]G^Q, 
and  Q,G^[Q,r]  =  QG^[r,P]  =  QG^rP,  as  QG-^P  =  0. 

[Q,r]G^Q  =    -   PrG^Q 

Q[G^r]Qr  =   -  |in  qg,  Z  Pij-G^Qr.  . 

(A59)  [Q  G^Q,r]r  =   QG-^rPr   -   PrG-^Qr    . 

Using   this   we   get 

<^WGr^GW7=  <WrGrGW^  +<'WGrPrGW> 

-  <WPrGrGW>  -  |^  <¥G  Z   p^jGr^jGW> 

Now   ^,    tp ,    W  and  G   are   even   functions,    and   r   is    odd. 
Therefore   the   second   term  vanishes    for    it   involves 
[^jMQfVip]    =   0.      The   third   term  Involves 

((t)W^)    =    (¥()), ^) 

=    ((E-H)(t),^) 

=    ((i),(E-H)^)    =    0    . 

In  effect,  [Q,,r]  has  been  set  equal  to  zero. 

<WGr'^G¥^  =  <WrGrG¥>-  |^   [<WG  Z   p^^.G   r^jW> 
-  <WG  2   p^jGQ[Gi,r^j]GW>     . 
Here  again  the   terms  with   [Q,r]    do   not    contribute. 
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<WGr^GW>  =  <WrG^rW>  -  ||^{-S  <WGr.'  .G^  p.  .GW> 

+  <WG  Z   p^jG   r^jW   +  ii^  <WG  Z   p.  jG^p.  jGW> 

h   M  Pi 

The  middle   terms   are  7^  —  Im  <WG  Z   p.  .G   r.  .W> 

(A60)     <;WGr^GW^  =   <WrG^rW>   +  -g  ^  Im  <WG  Z   p^  .G^r^  .   W> 

+  1    (|)2<WG  Z   p.jG2p^jGW>  . 

The  objective  of  commuting  r  through  has  been  obtained. 
We  can  express  matrix  elements  of  r.'  .  as  follows. 

'22,12,12  ,1 

r..  =  r..  +-n-r,  .  +-n-r..  -  r..  r,  .  -  r..  v„  .    +T^r,  .  r„. 


1 2    ^     2     ^   ' 
Z  r  .  =  ^  Z  rf  .  -  ^  Z  r.  .  r,  . 

_  3  .  9  ^2 
"22     • 

Now  since  X  is  symmetric  in  i  and  j,  we  have 

(a6i)         (^.r!^^)  =  |(^,r2/t:) 

Since  ^  depends  upon  relative  coordinates  only,  ^X=   0 
(A62)  p^-X  =  Z  p^X 
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13.   The  Relationship  Between  (f,Vf)  and  (f,V^f) 
for  Special  Choices  of  f. 

¥e  consider  the  function  f  =  f^e"^^-^  '     .      V^  Is 


defined  as 


V^  =  V^  2       ^ 


0  -  |x. I 

ay,  2,    4   2^ 
(a63)         (f,V°f)    =   3VQfQ    /   e       '^  ^  4Tr  X  dx  d_^ 

ay,  2^  4   2v 
(a64)                    (f,f)    =   fg    /   e       ^  5  dx  dy 


(a65)  (f,f)    =    (|^)^    (|)^/^   f 0  -  1    • 


(A66)  (f,V°f)   =  vg    (^)i/2 


ay,    2,  4   2v  2 

(a67)      -V^(f,Vf)  =  V;LVQ5f^e      °  ^  dx  d^,    g  =  ^ 

=  3V,V,  f^„(Z,J(a,g,-V2(4a^,-5/2 
=  3V,Vo   6-5/2  ^5/2    (I  ^  g,-5/2    . 

(a67')  -  V^(f,Vf)   =  JVj^Vq   (1+1  g)"5/2    . 

We  want 

(a68)  ^v^Vq  (1+1  g)-^/2  <  vg   y^^  • 

This    equation   Is 
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or 

(A69)  7^  a^  ^  V^  f^   . 

(a+6g)-^        "   ^ 

This  Inequality  is  satisfied  for  any  V-,  for  a 
sufficiently  small  or  large.   Suppose  we  maximize  the 
left  hand  side  with  respect  to  a: 

^ ^  =  [2a  -  -g!](a+6g)-^  . 

^^  (a+6g)^        ^^ 

This  is  zero  if 

2(a+6g)  =  3a 

or 

(A70)  a  =  12  g. 

This  value  of  a  yields  for  equation  (a69): 

lo^g 
We  have  an  equality  for 

(ATI)  V^  =  k7  —  -%  • 

V.  Is  the  same  as  equation  (43),  except  for  the  factor 

(A72)  ifM  -    ^-^5  . 

This  Indicates  that  our  choice  of  V-^  is  very  good. 
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14.   The  Alpha  Particle  Probl 


em. 


We  write  for  the  polarized  wave  function 


(AT3a) 


(A73b) 


^   =  ^  ^±j,Kl 


h   2,3,4  =^'l^ 


^(x^+a^y'^+b^z^) 

1 


J  could  be  J>}    but  It  turns  out  that  J  =  2  is  sufficient 
There  are  4j-l  variational  parameters. 
We  have  integrals  like 


^  ^^ij,k,r®^i'j',k',^'^  ' 


where  0  is  an  operator  symmetric  in  the  indices.   The 
indices,  being  dummy  variables,  can  be  relabeled  so  that 

(i'j',  k',i')  =  (12,3,4).   We  note  that  <t>ij^k,^  =   '^ji,k,^ 

4  ' 
This  leaves  us  with  -—^  =   12  terms .   We  write  them 

dlagrammatically  as  follows,  with  a  line  connecting 

the  first  two  symmetrical  indices: 


1     . 

1 

2, 

3, 

4 

2       . 

.        1 

2, 

4, 

3 

5      -^     . 

1 
1 

3, 
5, 

2, 

4 
2 

3   1,    2,    4 

4      . 

3  1,    4,    2 

5      • 

1 

4, 

4, 
4, 

2, 
3, 
1, 

3 

2 
2 

4   1,    2,    5 

6      . 

1 

4   1,    3,    2 

7      . 

5 

4   3,    1,    2 
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These  consist  of  all  12  permutations  of  1  j,  k,  H   with 
1  2  In  a  fixed  order.   Terms  1  and  2  have  weight  factor 
1,    and  terms  3  to  7  have  a  weight  factor  of  2.   We  have 
only  7  terms  to  evaluate  for  ((j),0(f)). 

We  have  the  following  relationships  between  the 
coordinates : 


(A74)  L2  ~  ^1  ^  ^ 


(A75)  £3  -  Zl  -  y  +  -^  2S 


To  get  rii-r,  ,  we  proceed  as  follows.   We  set 

r^i  -  r-,  =  ax  +  b  y  +  c_z  , 
where  a,  b,  c  are  numbers  to  be  determined.   We  have 

L}\-Li  =  ^L2   ~  ^^\   +  ^113  "  I (£1+1:2 ■*  "^  ^-4  ~  '|(-l+-2'^-3^ 

=  cr^  +  (b-^  |)r^  +  (a  -  |  -  ^)v^   -  I a+  |  +  |)r^ 
We  get 

(A76)  vj^    -Ill=z+-iy+|x 


(A77) 


x 

i.5  -  ^2  "  -^  "  2 


(A78)  Zi|-Z2=l+3l"l-' 

by  replacing  x  In  equation  (A76)  by  -x. 
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Zi|  "  -3  "  ^-4  "  1:2  ^"^^3 ".^2^ 


(A79)  zii.  -  Z3  =  z  - 1  y 


2   1-  'r.-R:^ 

—1  — 


^  -  ^  2  (r,-R) 


2 
¥e  want  to  evaluate  r   in  terms  of  x,  ^,    and  _z. 


64r2  =  [-x-(y+|)  -  (z+^  y  +|  x)]^ 


^        -      1   .2 


+  [x+(-2+  -^l  -2-31  +  2^^' 


+  tz+|+y-i2S-2+|y] 


+  [2+^y+|x  +  z+^-|  x+z-  |y]  2 


[  2x   +  I  y   +  z]  2   +   [  2x-  ^  y   -    z  ]  ^ 


+  [2  I  y-z]2  +  [3z]2 


We  get 


(A80)  r^  =    1    (x2   +^  y2   +1   z2) 

P 

We  now  want    the   quantity  >       r .  . . 

"Kj      ^^ 
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ZI  Zij  =  ^^  +  (Z+  I  2S)^  +  (z  +  3  Z  +  I  2S)^ 


KJ 


We  get 


+  (z  -  i  2S)^  +  U  +  3  y-  I  x)^+  (z-  I  y)^ 


i<j   ^-^  ^ 


or 


(a81)  T~  r?.  =  4r2  . 

2    2    2 
We  can  express  x  ,  y  ,  z   In  terms  of  r.,  i  =  l,3j 

as  follows 


X  =  ilg  "  £1 


Z  =  Z3  -  I  (£1+^2) 


-3 

4 


3  (Zl+£2^-3^  ' 


as  Vu    has  been  set  equal  to  -r., -rp-r  . 
We  get 

(A82a)  x^  =  r^  =  ^Zi  "1:2  ^  ^2 


(A82b)   y^  =  ^  ^1  +  I  £i-Z2+  i  ''2+''5"  -1-3  "  -^2*^3 


(A82c)   z^  =  ^  (r^  +  2r^-r2  +  rl+   r|+  2r^-r^+  Sr^-r^) 


If  we  set  A  ^IZ  £?  •  =  7  i^^+  4  y^  +  |  z^) .  we  have 
(A83)  7  =  2A  . 
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A  is  Laskar's  variational  parameter. 

Me  now  want  to  evaluate  the  kinetic  energy.   As 
shown  at  the  beginning  of  this  section,  we  need  only 


'  *12,5,t 


T  =  -  kV^  ,    k  =  -fi^/an.  Let  8.  ,  s  -^ — 


J 


+  52j(2x-a.y-  |  b .z)  .^S^ .(2a.y-  |b.z)  . 


+  a^.(2b.z).](j)^2,5,4: 


(A84)        !jl2.3^i.||Z  (4-H5a,+  |b,) 


_  k(2^)2  (8x^6a2y2,-i|b2z2) 


The  volume  element  of  our  integrals  is  either 
dx  d^  d£  or  dr^  dr^  dr  ,  with  -r^  =  r^   +  r ^  +  Z^- 

15,   Various  Integrals  for  the  Alpha  Particle. 

We  wish  to  evaluate  the  integral  1,    which  is 
defined  as 

r        2  2    2 

(A85)   I(d^,...,dg)  =  /  expL-d^r^-dgr^-rg-d^r^-di^r^ 

-  d^3_-r:^-dgr2-r^]  dr^  dr^  dr^ 
We  complete  the  square  as  follows.   We  note  that 
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p  Pi       p        p? 

a.x      +  p^x  =    (  ^  X  + ^    )      -        "^    . 

We   first  use   this   with  a-,    =   d.,  ,  6^    -   d^r^   +  d^r^    . 

X  1  ±  id — id  5 — ^ 

We  get   for   I 

_^/p    r  1 2      P-|  p  p 


p 
We  now   set  Qq  =   d^   -   do/'^d^,      p^  =  '^f>£.-^  "   ^p'^i-E^/^d^  . 

We  get 

2  P 

2       ,2      pS  o  df 


I  =    (a^Qp)    -^^     J    exp    [-r^    -r^  +  i^^^ ^3^^i|"  Tfd    ^  ^  ^Zi^Zp^Z^ 

2  2 

p 

Q/P   ^  P  d  d   d       -3/2 

(A86)    I  =   7r^/'^4^[(4d^d^-d2)(i^di^-  ^) -(2d^dg-d2d^)  f  2dg-  -f^)] 


We  now  want  the  integral 
(A87)   J  =  /  ^E;j_r^  +  ...  +  E^g-r^)  e   ^  ^      "^^  "^^Hl^IIp^Z^ 


2  ,      ,  _   .  ^   -dir^-...-d^2'i:3. 


J=  -ZlEi^l(d,,...,dg)  . 
1=1      1 

To  express  J  in  a  concise  form,  we  define  the 
following  quantities 

(A88a)  I  =  TT^/^  I^ 

(A88b)  ^1  =  ^  (WX-YZ)"^/^ 

The  result  for  J  is 


3  ^9/2  ^5  ^ 


(A89a)  J  =  3f  ^^'  II  JZ  E^  Fi 


1=1 


(A89b)   F^  =  ¥  d^/d^  +  X  4d^  -  Ydgd^/d^   -  Z  2dg 


(A89c)   F2  =  -  W  2d2/d^  +  Y  d^/d^  +  Z  d^ 


(A89d)   F  =  4  W 


(A89e)  Fi^  =  4  d-j_  X 


(A89f )   F^  =  -  2  d^X  +  Y  d2/d^  +  Z  dg 


(A89g)   Fg  =  -  2  Y  -  2  d-L 


(A89h)   ¥  =  4d-|_d   -  dg 


(A89i)  X  =  4di^  -  d^/d^ 


(A89j)  Y  =  2d^dg  -  dgd^ 
(A89k)   Z  =  2dg  -  dgd^/d-^ 
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